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Preface 


The work described in this report was performed under the cognizance of the 
Guidance and Control Division of the Jet Propulsion Laboratory. 

The purpose of the report is to document and disseminate information of prac- 
tical value to the engineer charged with the responsibility for developing an 
attitude-control system for a flexible space vehicle. It is essentially a treatise on 
the dynamics of flexible vehicles as viewed by the control systems engineer, and 
deals with control system design only to the extent of suggesting ways in which 
vehicle flexibility can be accommodated in both preliminary design and final 
evaluation of the attitude-control system. 

The major portion of the report is concerned with the derivation of equations 
of motion for the computer simulation of a controlled and nonrigid vehicle, and 
the development of coordinate transformations that facilitate simulation. Three 
basic approaches to this simulation are covered: discrete-coordinate methods 
(Section II), hybrid-coordinate methods, (Section III), and vehicle normal- 
coordinate methods (Section IV), 

The discrete-coordinate methods of Section II involve few restrictions or ap- 
proximations, and in some cases are as general as Newton’s laws for the simulation 
of the dynamic response of a collection of interconnected rigid bodies. The limi- 
tations of these methods stem primarily from the difficulty of creating the required 
mathematical model of a real vehicle without exceeding the practical limits im- 
posed on computation by considerations of budget, schedule, and computer 
capacity. 

The hybrid-coordinate methods of Section III receive the greatest emphasis 
in this report, both because they appear to be most useful and because they 
are the least familiar. These methods may be applied only when some portions 
of the vehicle (flexible appendages) undergo deformations that may reasonably 
be assumed to remain “small,” thereby permitting the transformation to modal 
coordinates for vehicle appendages. The key feature of this approach, as opposed 
to the discrete-coordinate method, is the possibility of truncating the matrix of 
modal coordinates. 

The vehicle normal-coordinate methods of Section IV involve transformations 
of all the kinematic coordinates of the simulation, and not merely the appendage 
deformation coordinates. These methods are accordingly more limited, and even 
when applicable, they may require more complex coordinate transformations than 
the hybrid-coordinate methods would involve. In the simplest cases, however, the 
vehicle normal-coordinate methods probably afford the most efficient simulation, 
since they then permit the most severe coordinate truncation. 

In addition to the material in this report for use in the simulation of a controlled 
flexible vehicle, there is material of value in preliminary control system design. 
Because of the modal coordinates employed in both the hybrid-coordinate methods 
and the vehicle normal-coordinate methods, the corresponding equations of motion 
permit simulations of varying dimension and complexity, depending on the degree 
of truncation. As shown in Section V, one can truncate so severely as to represent 
a very complex structural appendage by a single modal coordinate, and then 
manually calculate dynamic transfer functions for use in preliminary control 
system design. 
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Abstract 

The purpose of this report is twofold: (1) to survey the established analytic 
procedures for the simulation of controlled flexible space vehicles, and (2) to 
develop in detail methods that employ a combination of discrete and distributed 
(“modal”) coordinates, i e., the hybrid-coordinate methods. 

Analytic procedures are described in three categories: (1) discrete-coordinate 
methods, (2) hybrid-coordinate methods, and (3) vehicle normal-coordinate 
methods. Each of these approaches is described and analyzed for its advantages 
and disadvantages, and each is found to have an area of applicability. 

The hybrid-coordinate method combines the efficiency of the vehicle normal- 
coordinate method with the versatility of the discrete-coordinate method, and 
appears to have the widest range of practical application. 

The results in this report have practical utility in two arear: (1) complex digital 
computer simulation of flexible space vehicles of arbitrary configuration subject 
to realistic control laws, and (2) preliminary control system design based on 
transfer functions for linearized models of dynamics and control laws. 
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Dynamics and Control of Flexible Space Vehicles 


I. Introduction 

A. Motivation 

The development of an attitude-control system neces- 
sarily involves a dynamic simulation of the vehicle being 
controlled, but the accuracy required of that simulation 
may vary greatly from one application to another. As 
long as the attitude-control accuracy requirements are 
low, and the vehicle is relatively rigid, the “dynamics 
block” in a control system block diagram is easily gen- 
erated. Modem space vehicles are far from rigid, how- 
ever, and attitude-control accuracy requirements are 
increasingly stringent, particularly for optical observa- 
tions from space vehicles. Further improvements in con- 
trol system performance depend in many instances on 
improved simulation of the vehicle dynamics. The im- 
portance of this subject is reflected in the existence of 
die NASA Space Vehicle Design Criteria Monograph, 
Effects of Structural Flexibility on Spacecraft Control 
Systems (Ref. 1). This document includes relevant case 
histories and references, as well as design recommenda- 
tions. 

The incorporation of vehicle flexibility into control sys- 
tem design is a requirement neither new nor unique to 
space vehicle applications. Designers of control systems 


for missiles and aircraft have long ago found it necessary 
to consider vehicle flexibility in control system develop- 
ment. To some degree, the concepts developed in this 
earlier context can be applied intact to space applications 
(Refs. 2-10), but in many respects space vehicles are 
unique, and new concepts must be developed to deal 
with them. 

The environmental history of every artificial satellite 
or spacecraft is marked by a brief interval of vigorous 
acceleration and vibration during boost, followed by pro- 
longed functioning in a quiescent mode of operation 
characterized by extremely small loads and accelerations. 
Some missions may require a second period of violent 
acceleration for all or part of the system, but still there 
are two distinct and radically different dynamic environ- 
ments to be considered. The universal solution to the 
dilemma this poses for the structural designer has been 
the adoption of lightweight (and extremely flexible) de- 
ployable appendages. The resulting vehicle is relatively 
compact and rigid during the launch phase of its history, 
but after boost termination it emerges like a butterfly 
from its cocoon, extending antennas and booms and un- 
furling solar-cell arrays until the structure has undergone 
complete metamorphosis. 
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Although modern space vehicles can be found in a wide 
variety of configurations, certain dynamic features are 
sufficiently common to be described as characteristic of 
space vehicles, as opposed to missiles, aircraft, or land 
and sea vehicles. Most current space vehicles can be de- 
scribed as the combination of one or more essentially 
rigid bodies with one or more extremely flexible bodies. 
This natural separation of the structural subsystems of 
the vehicles into two distinct regimes is a consequence 
of the widespread use of deployable appendages, which 
are much less commonly found on surface or airborne 
vehicles. 

The structural subsystems of a spacecraft or satellite 
are often required to undergo substantial relative motions 
during mission performance, while large antennas, solar- 
cell arrays, instrument packages, or propulsion devices 
change their relative orientation. With the exceptions of 
variable-sweep and rotary-wing aircraft, most surface and 
air vehicles have dynamically significant moving parts 
only as semirigid rotors in the form of wheels, propellers, 
and propulsion subsystems. A rigid, symmetric, fixed-axis 
rotor is easily incorporated into a dynamic simulation of 
a vehicle, but the presence of a fully articulated flexible 
antenna on a space vehicle necessitates major changes 
in the formulation of its equations of motion, particularly 
when the relative motion of vehicle and antenna is sub- 
ject to closed-loop, nonlinear control. 

A different class of problem is introduced by the use 
of discrete damping devices in spacecraft subsystem 
vibration isolation or passive spin-stabilization nutation 
attenuation. 

All of these characteristic features of space vehicles 
present problems in dynamic simulation, and some of 
these problems are quite difficult to resolve by applica- 
tion of the methodology of aircraft and missile control 
system design. There muse be a comprehensive examina- 
tion of the question of dynamic simulation for attitude 
control of space vehicles, and new approaches must be 
developed for applications that are beyond the scope of 
the efficient utilization of traditional methods of analysis. 

B. Scope 

The present study is concerned with the development 
of methods for the dynamic simulation of flexible space- 
craft. The emphasis in this report is on analytic simula- 
tions, although it is recognized that analytic results must 
generally be confirmed or corrected by experimental mea- 
surements on models or prototype vehicles. 


Analytic methods treated here are restricted in their 
applicability to vehicles admitting of idealization as com- 
binations of rigid bodies, particles, continuous elastic 
bodies, and (in special cases) fluids. The equations of 
motion of a continuous mechanical system (with con- 
tinuous spatial variation of mass and flexibility) are par- 
tial differential equations. It is assumed in this report 
that any such equations descriptive of linearly elastic 
solid subsystems have been subjected to a “modal co- 
ordinate” transformation (as defined later in this section) 
and the resulting coordinates have been truncated to 
permit representation of system deformation with a finite 
number of modal coordinates. Thus continuous linearly 
elastic solids are admissible only in the sense that they 
can be represented either as a collection of interconnected 
rigid bodies or by a finite number of modal deformation 
coordinates. Fluids are admissible under three conditions 
only: (1) the fluid in a given container can be idealized 
as rigid, (2) the fluid can be represented by a finite 
number of modal coordinates, or (3) all aspects of the 
fluid dynamics can be ignored, except possibly the in- 
fluence of the fluid on energy dissipation. In consequence 
of the formal exclusion of continuous mechanical elements, 
the equations of motion are ordinary (and not partial) 
differential equations. 

Three distinct approaches to the analysis of flexible 
vehicles are treated in this report: (1) the discrete- 
coordinate formulation, (2) the hybrid-coordinate formu- 
lation, and (3) the vehicle normal-coordinate formulation. 
These methods are considered in turn in Sections II, III, 
and IV, with emphasis on the second method. A brief, 
qualitative description of each of the three general 
approaches is given in the following paragraphs. 

The vehicle normal-coordinate method is the traditional 
approach to the vibration analysis of elastic systems. This 
method is well documented in texts (Refs. 11-13), and 
space applications can be found in technical reports and 
journal papers (Refs. 2-10), so in this report it is described 
only insofar as this seems necessary for completeness of 
the study and for comparison with less familiar methods. 

In brief, this method consists of formulating equations 
of motion, whenever possible, as a system of independent 
(uncoupled) scalar second-order differential equations. 
For limited motions of a restricted class of physical sys- 
tems, systematic procedures can be generated for the 
determination of the transformation necessary to change 
from some arbitrarily selected coordinate system to the 
coordinate system corresponding to uncoupled scalar 
equations of motion. In general, the new uncoupled co- 
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ordinates (called normal coordinates ) do not correspond 
individually to the translations and rotations of discrete 
points or rigid bodies of the vehicle; instead, each coordi- 
nate is associated with a motion in which the entire vehicle 
participates. Coordinates that correspond to motion of 
more than one particle or rigid body of the system are 
called distributed or modal coordinates, as opposed to 
discrete coordinates. When all the dynamic equations of a 
system are uncoupled (as in the vehicle normal-coordinate 
method), the modal coordinates are called normal-mode 
coordinates for that system. Because the equations of 
motion are uncoupled, the vehicle can undergo motion 
in which only one of the scalar coordinates of the system 
participates. In this mode of motion, all points of the 
vehicle oscillate at the same frequency (the normal-mode 
frequency) and the vehicle undergoes periodic deforma- 
tions into the same deformed shape (the normal-mode 
shape). Each modal coordinate defines the amount of the 
response in the corresponding mode, suitably normalized. 
The independence of these normal-mode coordinates per- 
mits the independent calculation of their participation 
in the vehicle motion. This is the key feature of the normal- 
coordinate approach, because it permits the exercise of 
engineering judgment in determining which coordinates 
are so significant as to warrant retention, and which may 
be abandoned in coordinate truncation. 

Even this classical approach to vibration analysis has 
not been very widely used in spacecraft control system 
design and evaluation, so aspects of the method peculiar 
to this area of application are emphasized in this report. 
The assumptions underlying the vehicle normal-coordinate 
approach, and die theoretical and practical limitations on 
its use, are brought into focus. Appreciation of the restric- 
tions implicit in this approach is particularly important in 
space vehicle application, because the method breaks 
down when the system includes nonlinearities, rotors, 
discrete dampers, or articulated moving parts, and these 
are precisely the features just described as characteristic 
of space vehicles. Formulating equations of motion as 
first-order (state) equations eliminates some of these 
obstacles, but the practice of working with normal-mode 
coordinates for the entire vehicle is still quite restricted. 

All of these difficulties can be accommodated by avoid- 
ing modal coordinates entirely, relying upon a complete 
set of equations of a collection of interconnected rigid 
bodies considered to be representative of the vehicle. 
Because each rigid body of the system model, is discrete, 
and the coordinates employed in this approach are coordi- 
nates of position and/or attitude of the individual bodies, 


this approach is sometimes called the discrete-parameter 
or discrete-coordinate method. 

The discrete-parameter approach to space vehicle sim- 
ulation has received great emphasis since 1965, primarily 
because of its generality. A growing body of literature 
on this subject is becoming available (Refs. 14-22), and 
digital computer programs based on these and other 
formulations are becoming commonplace tools of analysis. 

The approach most frequently adopted (Refs. 14—21) 
involves direct application of the Newton-Euler equa- 
tions of translation and rotation to various subsets of 
bodies in the assembly (often to each of the bodies indi- 
vidually). Alternatively, Lagrange’s equations are some- 
times applied (Ref. 22). 

When application is restricted to a vehicle model com- 
posed of n rigid bodies interconnected at n — 1 point 
contacts so as to form a “tree” structure topologically, 
the Newton-Euler equations allow dramatic simplifica- 
tion. These restrictive assumptions are illustrated in 
Fig. 1, which shows that adjacent rigid bodies share at 
least one common point, and no closed loops are formed 
by any string of rigid bodies. The 12 bodies and 11 points 



Fig. 1. Topological tree of point-connected rigid bodies 
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of contact in the figure are arbitrarily labeled. The simu- 
lations described in Refs. 14-19 and 21 are limited in 
this manner (the “tree” of Ref. 21 is more like a “bush,” 
with f , trunk and limbs each consisting of a single rigid 
body). Under these conditions, for a complete kinematic 
description of the system, it is sufficient to know the 
location of one of its points (e.g., the mass center) and 
the orientation or attitude of each of the bodies. It is 
therefore possible to reduce the number of equations 
required for the dynamic simulation, eliminating all 
translational coordinates (except those of the system mass 
center, which may be determined independently in many 
cases). 

Euler’s equations of rotational motion of a rigid body 
(say, the ith body) are three scalar first-order differential 
equations in the variables which are measure 

numbers of the inertial angular velocity of that body 
for a body-fixed vector basis. These variables, collected 
for all n bodies, generally constitute the unknowns of 
the final differential equations, which are therefore first- 
order equations. Because a set of angular velocity mea- 
sure numbers docs not in itself provide a complete 
kinematic description of the system, an additional set of 
first-order equations must be included to permit the de- 
termination of attitude by angular velocity integration. 
These kinematic equations may be expressed in terms of 
any desired set of attitude parameters, e.g., direction 
cosines or a subset of direction cosines, Euler parameters, 
Euler-Rodrigues parameters (Ref. 23), Euler angles (3, 1, 
3 rotations), or Tait-Bryan angles (1, 2, 3 rotations). The 
last three alternatives involve only three parameters (the 
minimum number) for the attitude description of each 
body, but all three-parameter sets are plagued by iso- 
lated singularities that make numerical computations 
impossible for certain attitudes (without coordinate trans- 
formation). Most spacecraft-simulation programs employ 
a larger set of attitude parameters, accepting the incon- 
venience of working with a set of coordinates interre- 
lated by one or more constraint relationships. There 
appears to be no overwhelming preference between a 
subset of the direction cosines (usually six in number) 
and the four Euler parameters (or Cayley-Klein param- 
eters or elements of a quaternion). 

If one of the bodies of the set has more than one point 
in common with another body (e.g., the two are con- 
nected by a line hinge), this may be represented in the 
simulation by an additional constraint equation. Proce- 
dures exist for the systematic incorporation of such “joint 
constraints” into the dynamic equations, with consequent 
reduction in dimension of the matrix differential equa- 


tion of the system. Analysts disagree, however, on the 
advisability of this operation for computational efficiency 
(Refs. 14, 18, and 20). 

When Lagrange’s equations are written in terms of a 
set of independent generalized coordinates (Ref. 24, arti- 
cle 26), constraint equations are absorbed automatically, 
and the equations of motion are obtained immediately 
as a set of second-order differential equations of mini- 
mum number. The restriction to the use of independent 
generalized coordinates may, however, be a handicap, 
since this precludes the adoption of Euler parameters or 
direction cosines for attitude description. As a result, 
such a conventional Lagrangian formulation must employ 
a coordinate system not entirely free of singularities. In 
practice, most (perhaps all) programs based on Lagrange’s 
equations use attitude angles as coordinates (rather than 
Euler-Rodrigues parameters), and this selection intro- 
duces computationally cumbersome trigonometric func- 
tions into the calculations. This may be a greater handicap 
than the presence of singularities, particularly if digital 
simulation is anticipated. This approach is also restricted 
to systems with holonomic constraints, although this is 
perhaps not as severe a restriction in aerospace applica- 
tions as it would be for land vehicle applications that 
include rolling wheels. 

One may increase the scope and flexibility of the 
Lagrangian approach in at least two significant ways. 
With the introduction of Lagrange multipliers (Ref. 24, 
article 87), the restriction to a set of independent gen- 
eralized coordinates is relaxed, and one may adopt any 
kinematically complete set of coordinates, providing only 
that constraint equations in the form of equalities (not 
inequalities) exist in sufficient number to offset the coordi- 
nate redundancy. (Unfortunately, aerospace vehicle sub- 
system constraints such as damper and gimbal “stops” 
are in the excluded class of inequalities.) This approach 
extends the scope to certain nonholonomic systems, and 
it permits the use of redundant attitude parameters (e.g., 
Euler parameters). The price paid is the added dimen- 
sion of the problem. This method does not seem promis- 
ing in comparison with the Newton-Euler method, and 
appears not to have been pursued for general space 
vehicle simulation. However, applications to restricted 
prob^ms can be found in the aerospace journal literature 
(Ref. 25). 

The second relevant modification of the Lagrangian 
formulation is more fundamental. It is possible to avoid 
the reliance in Lagrange’s equations upon generalized 
coordinates, which must be sufficient to establish fully in 
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themselves the configuration of the system at any given 
time. One may alternatively write Lagrange’s equations 
in terms of quasi-coordinates (Ref. 24, article 30), which 
arc quantities whose differentials may be written as linear 
combinations of differentials of generalized coordinates 
and time. Since the angular velocity measure numbers 
Wj, o»!,. <».j of a given ith body qualify as derivatives of 
quasi-coordinates, this Lagrangian approach can produce 
results with some of the qualities of the Newton-Euler 
formulation. 

A recent adaptation by Kane and Wang (Refs. 26, 27) 
of the quasi-coordinate formulation seems well suited to 
complex spacecraft simulation. The method is applicable 
to certain nonholonomic systems and to systems with 
redundant (and constrained) coordinates; yet it auto- 
matically eliminates nonworking constraint forces and 
torques. Although the approach by Kane and Wang 
appears to combine certain advantages of the Newton- 
Euler method and the Lagrangian method, this path to 
complex system simulation has not yet been taken to the 
point of a multipurpose computer program for machine 
computation of the response of arbitrary discrete param- 
eter systems. 

Section II contains an introduction to the formulation 
of the equations of motion for discrete-parameter systems. 
The text describes without derivation the results of Refs. 
14 and 15, and notes some of the features of Refs. 16-20. 
The method of Kane and Wang is also briefly described. 
In view of the functional and structural beauty of the 
equations developed in these references, it seems unlikely 
that improvement can be realized by additional indepen- 
dent derivations, and the analyst confronting the problem 
of simulating complex discrete-parameter systems is 
advised to acquaint himself first with the referenced 
literature. 

Discrete-parameter simulations are not without serious 
disadvantages. Satisfactory simulation of real vehicles 
may require a great many rigid bodies in the model. The 
resulting differential equations are then of high dimen- 
sion, and their digital solution may be plagued by high- 
frequency responses that are of no interest to the engineer. 
There is no mechanism for truncating the matrix of 
coordinates retained in the determination of vehicle re- 
sponse, so as to eliminate these high-frequency responses 
and reduce the dimension of the problem. 

When the system configuration demands the generality 
of a discrete-parameter simulation, this alternative must 


be adopted. But when all or part of the vehicle allows 
the assumption of small, linearly elastic deformation, the 
most efficient simulation is that which combines discrete 
coordinates with distributed (modal) coordinates, retain- 
ing the generality of discrete coordinates where neces- 
sary, but securing the computational advantages of 
(truncated) modal coordinates where possible. The result 
is a hybrid-coordinate system that permits accurate 
simulation of complex modern space vehicles with a 
minimum number of coordinates, and with irrelevant 
high-frequency oscillations removed from the integration 
output by modal-coordinate truncation. Although this 
method is a natural combination of the other two ap- 
proaches, it has received very little emphasis in the tech- 
nical literature to date (Refs. 28, 29). For this reason, and 
because this method seems to combine the computational 
advantages of modal analysis with the needed gener- 
ality of the discrete-coordinate approach, the hybrid- 
coordinate method receives major emphasis in this report. 

The question of control system simulation per se is 
included in this report only insofar as this issue is coupled 
to the primary question of flexible vehicle simulation. 

II. Discrete-Coordinate Methods 

A. Aug men ted- Body Methods 

Hooker and Margulies (Ref. 14) and Roberson and 
Wittenberg (Refs. 15, 16) have observed that when a 
system of n point-connected rigid bodies is assembled in 
a topological tree (as in Fig. 1), certain inertia-like terms 
naturally appear in combination in the individual equa- 
tions of motion of each of the rigid bodies in the set. 
These combinations admit of physical interpretation as 
the inertia dyadies (or tensors or matrices) of abstractions 
called the augmented bodies. Briefly, the ith augmented 
body consists of the ith body of the set together with 
certain particles (point masses) attached to each of the 
joints of that body. The point mass attached to a given 
joint of the ith body equals the total mass of all of the 
connected bodies located “outboard” of the joint. For 
example, the augmented body 9 of the 12-body system 
shown in Fig. 1 is illustrated in Fig. 2. The mass center 
of the augmented body is called the connection banj- 
center (or simply the barycenter). The inertia dyadic of 
the augmented body with respect to the corresponding 
barycenter is the term that appears in the equations. 

A variety of approaches might be considered in formu- 
lating the Newton-Euler equations of motion of a system 
of rigid bodies such as illustrated in Fig. 1. Perhaps the 
most direct approach is to isolate each of the bodies in 
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Fig. 2. Augmented body 9 and its barycenter B a 


turn and write the translational and rotational vector 
equations for that body, recognizing that neighboring 
bodies contribute to the applied force and torque. This 
procedure leads to the appearance of augmented-body 
inertia dyadics, so it is the method used here. It should 
be recognized, however, that this is not the only possible 
approach; one might, for example, exclude the equations 
of motion of the nth body and instead write the transla- 
tional and rotational equations of motion of the total sys- 
tem, or of any connected group of bodies within the 
system. Other possibilities are explored in Section II-B. 

For an individual body, say, the ith body, the transla- 
tional and rotational vector equations of motion may be 
written as 

F‘ = m ; a' t 

{ (1) 

T' = H' = I' • cb‘ + co' X I* • co’ j 

where F' is the resultant force and T‘ the resultant torque 
applied to the ith body, m,- is the body mass, a ; the mass 
center inertial acceleration, H‘ the body angular momen- 
tum, I 1 the inertia dyadic, and co' the inertial angular 
velocity. The quantities T', H‘, and I’ are all referred 
to the body mass center. In the notation used in this 
report, dot ( • ) over a vector indicates time differentia- 
tion in an inertial reference frame. 

Among the constituents of the applied force F‘ there 
are “hinge forces,” i.e., forces applied to the ith body 
by contiguous bodies at the points of contact. These 
hinge forces generally contribute also to the torque T' 
about the ith body mass center, and typically this is the 
only coupling mechanism between the translational and 
rotational equations (1). 

It is evident by inspection of Fig. 1 that knowledge of 
the attitude of every body of a point-connected set is 
sufficient for the complete determination of the system 


configuration and attitude; it is not necessary to know 
mass-center position coordinates as well, since these must 
follow from the kinematic constraints imposed by the 
joints. Consequently it must be possible to combine the 
vector equations (1) for all n bodies in such a way as to 
obtain n vector equations of rotational motion that are 
internally complete, without any coupling terms to addi- 
tional translational equations. This is accomplished by 
first solving the translational equations (I) for the hinge 
forces in terms of mass-center accelerations, and then 
obtaining from- kinematics the mass-center accelerations 
in terms of geometry and rotational coordinates. 

The specific manipulations that accomplish this reduc- 
tion of equations (1) to a dynamically complete set of 
half as many equations can be found in Refs. 14-18. In the 
work of Hooker and Margulies (Ref. 14), the manipulation 
is performed at the level of explicit vector equations for 
the individual bodies, while in the work of Roberson 
and Wittenberg (Refs. 15, 16), this task is deferred to a 
stage of the derivation in which the rotational equations 
of motion are matrix equations of dimension 3n. In both 
developments, there emerge terms (dyadics in Ref. 14 
and matrices in Refs. 15, 16) that have the indicated inter- 
pretation as inertia dyadics or matrices of augmented 
bodies. The vector-dyadic result, for example, is given 
for the Ath body in the notation of Ref. 14 by 

<D^ • cb* V co x x O*' c^a = 

Ta + 2 T". + Da X Fa + 2 Da m X F M 

jcA /i * A 

+ Qll 2 Pi/i X [cb M X D m ,\ + co M X (oj/x X Dfu)] (2) 

fi * A \ 

with the following symbol definitions: 

O* — abe inertia dyadic of the Ath augmented 
body, referred to the corresponding bary- 
center 

coxandw fl = i. ... inertial angular velocities of bodies A 
and /X., respectively 

Ta = that portion of the resultant torque ap- 
plied to body A obtained by excluding 
forces and torques applied at joints 

T",- = the “hinge torque” applied at joint j of 
body A 

/a = the set of numeric labels for the joints on 
body A (e.g., from Figs. 1 and 2, /„ includes 
7, 8, and 9) 
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Da = the position vector from the baryeenter 
B\ to the mass center of body A (e.g., see 
Fig. 2) 

F,\. and = the forces applied to bodies A and p, re- 
spectively, excluding forces applied at 
joints 

Da(.i = the position vector from the baryeenter B\ 
to the joint of body A that leads to body p 
(even if body p is not directly connected 
to body A, but instead is part of a chain of 
bodies connected to body A) 

Q 11 = the total system mass 

2 = the sum over values of / in the set /a. 

j eJ\ 

The left side of Eq. (2) has exactly the form of the 
vector-dyadic representation of Euler’s equations, except 
for the substitution of the augmented-body inertia-dyadic 
about the baryeenter (O*) for the body inertia-dyadic 
about the mass center (l A ). The right side of Eq. (2) in- 
cludes the body torque T v and the relevant hinge torgues 
T(1 , as would be the case for Euler’s equations, but the 
torque contribution of hinge forces takes a surprising 
form. In every case, torques are calculated with respect 
to the baryeenter, and the appropriate force turns out to 
be the sum of the body forces F ;i and certain “inertial 
forces” that may be attributed to angular accelerations 
and centripetal accelerations. 

In view of the complexity of the n-body system under 
examination, it is astonishing that the equations of motion 
should be so simple in structure and amenable to physical 
interpretation. Simplification is even more dramatic when 
there is a coincidence of baryeenter Ba and a joint, since 
then the vector Daju is zero for all bodies connected to 
body A by that joint, and substantial decoupling of the 
equations results. 

In general, of course, there are required n vector-dyadic 
equations such as Eq. (2) to complete the dynamic de- 
scription, and in addition some specification of the hinge 
torques TJj is required. In most applications, some of 
the bodies in the model are connected by line hinges, 
so that a line is common to two bodies of the system. 
The hinge torques transverse to the hinge axis then 
become constraint torques, and each line hinge provides 
two scalar constraint equations that preclude relative ro- 
tations except about the hinge axis. Similarly, a two- 
gimbal joint provides one constraint equation. If there 


are v constraint equations and n rigid bodies in the sys- 
tem, then there are available 3 n + v scalar equations to 
be solved for the 3 n unknown angular velocity measure 
numbers plus the v constraint torques. (An additional 
equation is of course required for the determination of 
the translation of the system mass center, but this equa- 
tion is generally uncoupled with the attitude equations 
and is ignored here.) 

The most direct approach to the determination of 
dynamic response is to combine all the equations of dy- 
namics and kinematics with the constraint equations as 
a single-matrix first-order differential equation of dimen- 
sion 6n + v and proceed with numerical integration. 

The possibility of eliminating the v constraint torques, 
thus reducing the dimension of the problem to 6n, is 
discussed briefly in Ref. 14, and related questions are 
explored er tensively in Refs. 16, 18, 19, 20, 26, and 27. 
Some of the methods that serve this purpose are discussed 
in Sections II-B and II-C. 

The task of assembling the n vector-dyadic rotational 
equations into a single matrix equation is made quite 
awkward by the summations over limited sets of joints 
and bodies that appear in Eq. (2). Roberson and 
Wittenberg (Refs. 15, 16), in a derivation that developed 
in parallel to that of Hooker and Margulies (Ref, 14), 
utilize graph theory to construct an elegant formalism 
for the systematic assemblage of the individual vector- 
dyadic equations into the single 3n-dimensional matrix 
equation of motion of the total system. The equations of 
Ref. 15 are programmed (with some modification) in 
Ref. 17, and those of Ref. 14 provide the basis for the 
program developed in Ref. 18. 

Whether the matrix equation is developed with the 
beautiful formalism devised by Roberson and Wittenberg 
or assembled in more pedestrian fashion (as computa- 
tional efficiency may suggest), the final structure of the 
equation must be as follows: 

ft o> = m + j ? ( 3 ) 

where T’ absorbs all terms involving constraint torques, 
0) embraces all external forces and torques as well 
as all terms involving inertial angular velocity terms 
o)' a (a = 1, 2, 3 , i = 1 , • • * , «), and ft is the coefficient 
matrix of all differentiated terms, which are assembled 
in the column matrix «. Ignoring for the moment the 
question of how constraint equations are to be used to 
accommodate unknowns in -i?, one can readily see that 
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it is the nature of (ft that will determine computational 
difficulty. If (ft were a nonsingular diagonal matrix, a 
trivial inversion would permit Eq. (3) to be written as 

+ ( 4 ) 

where superscript —1 denotes matrix inverse. 

Unfortunately, a glance at Eq. (2) indicates that iner- 
tial coupling terms involving must exist unless D,\ M = 0 
for all /a. The same term D,\ M X c also assures the time- 
varying character of (ft, since the cross product depends 
on the changing relative attitude of bodies A and /e This 
constitutes a major obstacle to numerical integration, 
since if Eq. (4) is used it becomes necessary to invert 
the 3 n by 3 n matrix (y-\ at each step of the integration. 
In practice, it may prove more convenient to retain the 
equation in the form of Eq. (3), applying Gaussian elimi- 
nation rather than matrix inversion at each integration 
step. Even with this expedient, the algebraic process of 
finding & is apt to consume most of the computer time 
in numerical integration. Since time required for this 
process increases roughly with the third power of the 
dimension of the matrix (ft, abundant motivation exists 
for working with equations of minimum dimension. 

B. Nested-Body Methods 

In writing n vector equations of rotational motion for 
a set of n point-connected rigid bodies, one must make 
a choice of the n material subsystems to be isolated for 
dynamic analysis; it is not obviously advantageous to 
isolate each of the rigid bodies individually, as is gen- 
erally the practice in Refs. 14-18. (In Refs. 15-18 there 
is an explicit departure from this practice when a pair 
of rigid bodies consists of a rigid body containing a rigid, 
axisymmetric rotor on a fixed axis.) 

Velman (Ref. 20) and Russell (Ref. 19) elect instead 
to write vector equations of motion in turn for n different 
subsets of bodies, including a final set of rotational and 
translational equations for the composite vehicle. For 
example, if the vehicle were modeled by 12 point- 
connected rigid bodies as shown in Fig, I, both Velman 
and Russell might write equations of motion first for 
body 12, then for the subset including bodies 11 and 12, 
and then for 10, 11, 12. Next, body 8 might be isolated 
and its equations written, and these followed by equa- 
tions for the subset including bodies 8-12. One might 
then begin anew at the end of another chain of bodies 
(either body 1 or body 7), and progress inward as far 
as possible without ever including more than one addi- 


tional body in any one step, and without ever considering 
a subset with more than one joint connecting it with the 
excluded bodies of the system. For this example, one 
might follow the indicated pattern until each of the three 
chains attached to body 4 has been considered, without 
ever including body 4 itself. (This choice of body 4 is 
arbitrary, since other subset selections within the pre- 
scribed pattern can converge as well on any other body 
of the system.) Finally, the equations of motion of the 
composite vehicle are recorded, to complete the dynamic 
simulation. 

The idea of isolating in sequence such subsets of rigid 
bodies in an n-body system (called nested bodies by 
Velman) seems to have both advantages and disadvan- 
tages. The concepts of the augmented body and the 
connection banjeenter are helpful aids to physical inter- 
pretation (if not to numerical computation), and the 
terms with which these phrases are identified in the 
equations of Roberson- Wittenberg (Refs. 15, 16) and 
IIooker-Margulies (Ref. 14) do not appear when the 
nested-body approach is taken. In compensation, how- 
ever, it would appear that the nested-body derivation 
facilitates the elimination of internal constraint forces 
and torques. 

Velman (Ref. 20) derives his equations in a manner 
that leads quite naturally to the use of attitude variables 
that establish the orientation of each body of the system 
relative to an adjacent body (except of course for one 
reference body with inertial attitude parameters). Because 
geometrical constraints generally restrict relative motions, 
it is perhaps easier to identify which attitude variables 
are constrained when relative motion coordinates are 
employed. Nonetheless, there is no attempt to eliminate 
kinematically constrained variables in the derivation of 
Velman’s dynamical equations, so his equations have at 
the outset the same basic structure and dimension as those 
of Roberson and Wittenberg or Hooker and Margulies, 
namely, 

(ftu> = m + j? (5) 

These symbols differ in interpretation among the several 
authors, and, in particular, Velman’s w includes some rela- 
tive velocities. It should be noted also that Velman’s sim- 
ulation specifically includes a linearly oscillating particle, 
so the method is not restricted to a point-connected set 
of rigid bodies. In Eq. (5), then, one of the scalars in the 
matrix « is the translational velocity of a particle relative 
to its housing body. As in Refs. 14-16, the dimension of 
matrix <i> is 3 n for an n-body system as long as each body 
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is either a particle or a rigid body forming part of a point- 
connected set. Again, A} and 0} are matrices that depend 
upon the kinematic unknowns of the problem, and -C 
is the matrix of unknown constraint torques (and con- 
straint forces for the oscillating particle). The most sig- 
nificant feature of Velman’s work is the procedure he 
devises for the elimination from these equations of the 
unknown constraint forces and torques. 


9 


0 

0 

0 

0 

0 

0 


0 0 
0 0 
0 0 
0 0 
0 0 
0 0 


0 0 0 “ 
0 0 0 
0 0 0 
10 0 
0 10 
0 0 1 


(7) 


Velman notes that the effect of the constraints is usu- 
ally to confine the solution of Eq. (5) to some linear 
manifold in the 3/i-dimensional space of u>. If, for exam- 
ple, a line hinge connects two rigid bodies of the system, 
and three of the scalars in « are the relative rotation rates 
of these bodies about the hinge axis and two transverse 
axes, the effect of the constraints is to confine the solu- 
tion to that subspace of the o> space excluding the two 
transverse-axis rotation rates, which are constrained by 
the hinge to be zero. The subspace to which the solution 
is confined is then of dimension 3 n — 2. This fact is 
unchanged even when the choice of variables in w does 
not specifically include these three relative rates about 
and transverse to the hinge; this choice may simplify the 
explicit specification of the solution subspace, but it intro- 
duces no conceptual change in the ary ament. Similarly, 
the line of argument remains intact when the solution 
in a subspace of «> is constrained to be not zero, but a 
specified function of time. If, for the previous example, 
the two bodies connected by a line hinge undergo a pre- 
scribed relative rotation about that hinge (as for an ideal- 
ized scanning antenna), then the solution is a prescribed 
function of time in a three-dimensional subspace of the 
o> space. 


The preceding remarks amount to the observation that 
the constraints determine the projection of the motion 
on a particular linear manifold in the w space. Velman 
introduces the perpendicular projection operator 9 (a 3n 
by 3/i matrix called K in Ref. 20) and the analytic expres- 
sion 

<? o> - (0 (6) 

to represent the specified motion /<. (f), which is the con- 
strained partial solution. Because the variables in w are 
generally the relative motions, the matrix ( ? is typically 
a diagonal matrix, null except for ones on the diagonal 
corresponding to constrained coordinates. For a two- 
body system consisting of a primary body with inertial 
rates i.) b w 2 , and a second body attached on a line 
hinge with prescribed relative rotation rates w, = 0, 
CD.-, = 0, wr, — 0 (t), the ( P matrix would have the structure 


and the kinematic constraint equation would be 
9 w — [0 0 0 <i).j <!>-, wg]^ 

= [00000 e ( t)] T ( 8 ) 

where superscript T denotes the matrix transpose. 

Rather than augment the 3n scalar dynamic equations 
in Eq. (5) by the v scalar constraint equations that con- 
stitute the nontrivial part of Eq. (6), Velman (Ref. 20) 
uses the latter to eliminate the constraint torques 1’ 
from the former. To accomplish this, he notes that the 
vector of constraint torques (and forces) -C lies in the 
same manifold of the co space in which motion is pre- 
scribed. Thus if -7 is defined as the supplementary per- 
pendicular projection operator for 9, so that 9 + 9 is 
the identity matrix E, 


' 7 + 9 = E 

then 7 operating on 1’ is given by 

9J? = 0 


(9) 


( 10 ) 


For the simple example of the two hinged bodies, Eq. (10) 
appears in expansion as 


"1 0 0 0 0 0" 


" 0 ' 


'"0“ 

0 1 0 0 0 0 


0 


r 

0 0 1 0 0 0 


0 


0 

0 0 0 0 0 0 


TV 


0 

0 0 0 0 0 0 


,XV. 


0 

_0 0 0 0 0 0_ 


_TV_ 


_0_ 


( 11 ) 


Now one can return to Eq. (5) and rewrite it as 


ffl(9 4- 7) oj = 0} +J? (12) 
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which upon premultiplication by 0 and substitution of 
Eqs. (6) and (10) becomes 

0919 « = 0m - 991 IX (: t ) (13) 

Velman notes that the matrix 99l9 is singular, so 
Eq. (13) is not yet ready for machine computation. Using 
the idempotency of the perpendicular projection oper- 
ators and their supplementary character, one may write 

99 = 9, 09 = 9, 09 = 0 (14) 

and then rewrite Eq. (13) in the form 

(9019 + 0)9* = 9rc] - 9 91 ,, (t) (15) 

The matrix (99\9 + 9) can be inverted to obtain 

9 * = {9919 + 0)~' 9[m-m p- (t) ] (16) 

which when added to Eq. (6) yields 

o. = {9919 + 9)-' 9 [re -9lp(t)]+p (t) (17) 

since 9 + 9 = E from Eq. (9). This is the form of the 
equation used by Velman. The dimension of the unknown 
matrix m is still 3 n, but the constraint torques have been 
eliminated. 

As noted earlier, the fact that Velman’s variables in w 
are generally relative motions simplifies the physical 
interpretation of the constraint equation (6), and often 
results in constant perpendicular projection operators 0 
and 9. Still, the coefficient matrix 91 depends on the 
equation variables and accordingly varies with time, so 
the 3 n by 3n matrix {9919 + 9) must be inverted at 
each integration step for the dynamic equation (17). Since 
the system has only 3 n — v degrees of freedom, one might 
hope to reduce the dimension of the variable to in Eq. (17), 
and in the process reduce the size of the matrix requiring 
inversion at each integration step. 

Fleischer (Ref. IS) has developed a procedure for 
accomplishing the indicated reduction of the dimension 
of the problem to 3 n — v. If the variable u> is partitioned 
as w = [ iof * C ] T , with subscripts f and c identifying free 
and constrained coordinates, then one may attempt to 
use the constraint equations to find a> ( . in terms of *(, and 
vhen retain only that portion of Eq. (17) involving £>/■ on 
the left side (i.e., truncate the matrix w to <h r ). Actually, 
Fleischer pursues a somewhat different path from the 
basic equations (5) and (6) to the final equations of mo- 


tion, and in the process he introduces new operators to 
facilitate the reduction of dimension of the final equa- 
tions. Because the appropriate modifications are fully 
treated in Ref. 18, they are not included here. 

Fleischer’s final equations have the form 

9lf «V — 9 >[ (18) 

in which 9lf and are time-varying and of dimension 
(3» — v) by (3n — v). Computation requires the inversion 
of 01 f at each integration step (or the equivalent Gaussian 
elimination process), and since this inversion is the most 
time-consuming part of a major simulation program, this 
reduction of dimension would appear to be a substantial 
improvement. Yet there are additional computer opera- 
tions involved in obtaining and solving the reduced 
equations, and differences in the detailed structure of 
the equations may impede computation in the reduced 
case, so perhaps it would be useful to discuss the rela- 
tive computational advantages of the several forms of 
the discrete-parameter system equations. For example, 
buried within the matrices 9l f and ti\ t in Eq. (18) there 
is the inverse of a matrix of dimension v. When <.> com- 
prises relative motion variables (as in Velman’s deriva- 
tion), the matrix to be inverted is generally constant, 
requiring only one inversion operation. But when w com- 
prises inertial rates of the various bodies, this matrix 
depends on time, and the matrix inverse must be obtained 
with each integration step. In Fleischer’s work (Ref. 18), 
the basic dynamic equations of Hooker and Margulies 
(Ref. 14) are combined with a modified version of 
Velman’s constraint elimination procedure (Ref. 20), so 
inertial rates are the variables and the indicated repeated 
inversion is required. At each step of the integration, 
Fleischer must invert a v by v matrix and then invert (or 
apply Gaussian elimination to) a matrix of dimension 
(3n — v) by (on — v), whereas Velman must invert (or 
apply Gaussian elimination to) one matrix of dimension 
3 n by 3)). As noted, Hooker and Margulies appear to 
favor working directly with the dynamic equations as 
augmented by the constraint equations, accepting system 
equations of dimension (3 n + v) by (3)7- + v). It seems 
probable that selection of an optimum approach will vary 
from one application to the next, but that for very com- 
plex systems the advantage would lie with the approach 
yielding equations of the lowest dimension. 

In the derivations of Refs. 14-18 and 20, the approach 
has been first to write dynamic equations in which con- 
straint forces and torques appear, and then (in some 
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cases) to use constraint equations to eliminate from the 
dynamic equations the constraint forces and torques and 
the redundant kinematic variables. This elimination is 
not accomplished by the analyst before computation, but 
in general is accomplished by the computer during 
numerical computation. 

Other approaches have been devised that avoid con- 
straint forces and torques from the outset, and thereby 
preclude the necessity of using computer time to elimi- 
nate constraint forces and torques and reduce the num- 
ber of scalar equations to the number of degrees of 
freedom. In the context of the Newton-Euler formula- 
tion, the objective of avoiding constraint torques for a 
point-connected n-body system in a topological tree has 
been pursued successfully by Russell (Ref. 19). Kane and 
Wang (Refs. 26, 27) accomplish this objective even more 
generally, using a method described in Section II-C. 
Although neither Russell nor Kane and Wang support 
the concept of developing a single multipurpose com- 
puter program suitable for a wide range of vehicles, 
Russell docs provide in Ref. 19 an explicit procedure for 
constructing equations for simulation of an n-body sys- 
tem restricted as above. 

Russell adopts the nested-body concept advanced by 
Velman and described in Section II-A, where the 12-body 
system of Fig. 1 was discussed as an example. Recall that 
for this example one may first write the vector rotational 
equations of motion for body 12, then for subsystem 11 
plus 12, etc. As pursued by Velman, this path does 
involve constraint torques. 


(For this reason Russell calls his method the momentum 
approach.) This distinction is central to the computational 
question, but not critical to the avoidance of constraint 
torques. 

Having obtained a single scalar equation for the com- 
ponent of H 1 - along the hinge axis, one next proceeds to 
consider the subsystem of bodies 11 and 12, writing trans- 
lational and rotational equations as previously. Now the 
analyst must inspect the single joint between this subsys- 
tem and the rest of the vehicle, namely, the joint connect- 
ing bodies 11 and 10. If this is a 3-degree-of-freedom 
joint, the entire vector rotational equation of the sub- 
system 11 plus 12 is retained, but if this is a rotationally 
constrained joint, only those components of the vector 
equation that do not introduce constraint torques are 
retained. Whatever the nature of the joint, the number of 
scalar equations added corresponds to the number of 
degrees of freedom added in considering a new subsystem. 
Proceeding in this way, following the pattern of selecting 
subsystems previously described, one can systematically 
accumulate as many dynamic equations as there are inde- 
pendent unknowns without introducing constraint torques 
at all. A simple accounting procedure permits the distri- 
bution of the angular momentum of any subsystem among 
its constituent bodies, and from these individual inertial 
angular momenta the inertial angular velocities follow 
from Eq. ( 19 ) . The dynamic equations then combine with 
kinematic equations just as they do for the alternative 
Newton-Euler formulations (Refs. 14-18, 20). 

C. Generalized-Force Methods 


Following Russell’s approach, the analyst must observe 
that body 12 is connected to body 11 by a line hinge, and 
then write for body 12 only that component of the rota- 
tional equation T 12 = H 12 paralleling the hinge axis. This 
avoids the constraint torques (assuming unspecified rota- 
tion about this axis), and it yields one scalar equation 
corresponding to the single degree of freedom of body 12 
relative to body 11. Because the torque T 12 is measured 
relative to the mass center of body 12, the interaction 
force applied to body 12 by body 11 contributes to T 12 , 
and this unknown constraint force must be eliminated 
by use of the translational equation F 12 = m vi a 12 (Ref. 19). 
The method advanced by Russell retains the components 
of the inertial angular momenta of the individual bodies 
as the unknowns, obtaining angular velocity of the ith 
body as necessary from the matrix relationship 

a, 1 ' = (V)- 1 W (19) 


Among the many methods employing generalized forces, 
most are within the framework of Lagrange’s equation 


cl f ?L\ cL_ 
clt \?q,) 


■ ■ ,n f 


( 20 ) 


where q u • • • ,q„ f are a complete and independent 
set of generalized coordinates, the Lagrangian L is the 
difference in kinetic and potential energy, and the gen- 
eralized force Qi is defined in terms of applied forces 
F 1 , • • • ,F ,S ' and their inertial position vectors r 1 , • • • ,r ,s ' by 


Qi 


/.I rf/i 


i -!,•••, n f 


(21) 


Lagrange’s equation in the form of Eq. (20) is too familiar 
to warrant review and too limited to warrant adoption as 
the multipurpose procedure for discrete-parameter sys- 
tem simulation. As noted in the Introduction, restriction 
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to a complete and independent set of generalized coordi- 
nates q u ■ ■ ■ , q„ { is inhibiting in that it precludes the use 
of redundant, singularity-free attitude variables and ex- 
cludes nonholonomic systems. (A holonomic system can 
be fully described by a set of coordinates q !,•••, q x , which 
are related by v equalities of the form f t (q u • • • , q x , t) — 0, 
i — 1, • • • , v, thus permitting, at least theoretically, the 
determination of v coordinates in terms of the remaining 
set of iif — N — v coordinates.) 

The generalized forces Q-, are attractive, however, in 
that the definition in Eq. (21) serves to eliminate from Q-, 
any nonworking constraint forces, and thus to eliminate 
these unknown and unwanted quantities from the equa- 
tions of motion. 

Kane and Wang (Refs. 26, 27) have devised an approach 
that retains generalized forces with this desirable prop- 
erty, but permits application to certain nonholonomic 
systems and to systems with a wider range of acceptable 
coordinates. 

Consider a dynamic system described by a complete 
but redundant set of coordinates q u • ■ • , q x subject to 
the v- constraint equations typified by 

2 A }i qt + B } = 0, j = ( 22 ) 

< -- 1 

where Ay and B; are functions of q u • • • , q x and t. If 
Eq. (22) cannot be integrated to obtain constraint equa- 
tions in the form of equalities not involving generalized 
velocities, this is typical of a class of nonholonomic 
systems. 


In application to a system of n particles and rigid 
bodies, one may without restriction write the Newton- 
Euler equations for the ith body of the system in the form 


where 

F>' + F'* 
T‘ + T i!S 

II II 

o o 

II 

i— > 

■ • ,n 

(23) 

(24) 

and 


F 1 * = — nija 1 


(25) 



T- - — H‘ 


(26) 


with T', a', and H' referenced to the mass center for the 
rigid bodies. In application to a rigid body, F‘ is the 
resultant applied force. Because the constraints imposed 


in Eq. (21) imply the presence of constraint forces in F' 
and constraint torques in T', these vector equations of 
motion when applied directly do not have the constraint 
elimination feature of the Lagrange Eq. (20). The objec- 
tive of Kane and Wang is to modify Eqs. (23) and (24) so 
as to secure this advantage without accepting the noted 
limitations inherent in Lagrange’s equation in the form 
of Eq. (20). 

To apply the Kane- Wang method, one must first gen- 
erate the necessary kinematic quantities and record the 
inertial forces F'* and inertial torques T'* for all bodies 
of the system. This process includes the derivation of the 
inertial linear and angular velocities v' and £»>' of all 
bodies of the system, and their expression in terms of the 
constrained generalized coordinates q u • • • , q x and their 
first derivatives. It is always possible to use the v constraint 
equations given hyEq. (22) to write all velocities v' and co 1 
in terms of N — v (say, the first N — v) generalized veloci- 
ties q u • • • ,qx-v, retaining in general all N coordinates 
q i, • • • ,q x in these expressions. 

The next step is the selection (by inspection of the 
expressions for v' and <o') of N — v variables u u ' ' ' , Ux-v 
so that (1) each tq may be written as a linear combination 
of the IV — v generalized velocities (/,,••• ,q x - v plus a 
residual term free of generalized velocities, i.e., 

m = S’ Uu qj + U h i = 1, • • • , N - v (27) 

/= i 

(where Uy and 17; depend on q„ • • • , q x and t), and 
(2) each velocity v' and angular velocity W may simi- 
larly be written as a linear combination of the variables 
«i, • • ' , M.v-v, i-e., 

v 1 ' = 2 vj « j + vj (28) 

and 

•V - V 

co 1 ' = 2 00 j Uj + “i (29) 

This step clearly involves some judgment on the part of 
the analyst, but in specific applications the choice of vari- 
ables iq is not difficult. It is always possible to choose 
u\ — qi, i = 1, ■ • • , N — v, but this is rarely the most 
attractive choice. More often a judicious choice of angu- 
lar velocity measure numbers for iq, • • • , ttx-v is useful. 
Examples in Refs. 26 and 27 are helpful in establishing a 
rationale for this selection. 
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From this point, the application of the method is quite 
routine. The vectors vj, coj, vj, • • • coj, • • • , coj. t , 

(/ -- 1, • • • , n) are defined by the selection of variables 
u i, • • ■ , 11 . x-v in Eqs. (28) and (29). These vectors are 
now used to define the “generalized active forces” 

K; - i (v) * F' + coj • T'), / = !,••■ - v 

(30) 

and the “generalized inertia forces” 

K] = S (vj • F' - * + coj • T'*), j = 1, • • • ,N — v 

(31) 

Kane’s theorem (proven in Ref. 26) has the form 

Kj 4 K* = 0, / — 1, ' ' ' ,N — v (32) 

Kane shows that the basic first-order dynamic equa- 
tion (32) does not involve the unknown forces and torques 
required to maintain the constraints of Eq. (22). Thus the 
combination of the N — r dynamic equations of Eq. (32), 
the v constraint equations of Eq. (22), and the N — v kine- 
matic equations of Eq. (27) constitutes a complete set of 
2 N — v first-order differential equations. This is the mini- 
mum number of equations consistent with the selection 
of a system of N generalized coordinates subject to v con- 
straints. In application to the point-connected sets of n 
rigid bodies of Section II-B, N is 3n, so Kane’s equations 
are the same in number as those obtained by Russell’s 
approach (Ref. 19), and the same as those actually inte- 
grated in Fleischer’s program (Ref. 18). If the constraint 
equations are holonomic (so Eq. 22 is integrable), it is pos- 
sible to replace v of the first-order differential equations 
by algebraic equations. When these can be solved ex- 
plicitly for v coordinates in terms of the N — v remain- 
ing coordinates, the results can be substituted into the 
dynamic equation (32) and the kinematic equation (27) 
to obtain a minimum set of 2 (A 7 — v) first-order differ- 
ential equations. 

D. Summary 

The objective here lias been to outline the several 
approaches to the formulation of equations of motion of 
discrete-parameter systems in sufficient detail to develop 
some appreciation of the differences in points of view, 
and to stimulate interest in deeper inquiry. For a compre- 
hensive treatment of the subject, the reader is referred to 
the cited works. 


It is impossible to make a definitive judgment of the 
relative merits of the many available approaches. It is 
difficult to weigh computer time against analyst time, 
particularly in the face of absolute time constraints and 
limitations of personnel assigned to a given task at a given 
time. Even with a given “cost function,” the optimum 
approach varies from one dynamic system to the next. 

Within the framework of fully programmed multi- 
purpose digital computer programs, which require the 
minimum of analyst time and thought, the work of 
Fleischer (Ref. 18) is perhaps most useful, since it com- 
bines many of the attractive features of the earlier works 
of Hooker and Margulies (Ref. 14), Roberson and Witten- 
berg (Refs. 15, 16), and Velman (Ref. 20). It is possible 
that Fleischer’s program would run faster if relative 
motion coordinates were employed, but this would intro- 
duce the judgment of the analyst into the simulation 
process. 

There may be some question concerning the desir- 
ability in some organizations of developing a multipur- 
pose program. Russell (Ref. 19) and Kane and Wang 
(Refs. 26 and 27) offer procedures that seem to lead to 
more efficient simulations than are available with the 
more highly formalized methods cited above, but they 
require individual programming of each problem by a 
capable analyst, which may be an unacceptable constraint. 

III. Hybrid-Coordinate Method 
A. Vehicle Mathematical Model 

The concept of the hybrid-coordinate method depends 
upon the possibility of separating a given vehicle into a 
number of idealized structural subsystems, each of which 
may be classified either as a flexible appendage or as a 
rigid body or particle. A flexible appendage is by present 
definition a linearly elastic structure (developing elastic 
stresses in proportion to strains) for which “small” defor- 
mations may be anticipated (so elastic stresses remain 
proportional to deformations). In most cases, internal 
energy dissipation may be represented by modal damping 
(Ref. 30), although discrete dampers can be accommo- 
dated. The usual definitions of rigid body and particle 
are retained. 

In its present stage of development, this method is 
formally restricted by the requirement that each flexible 
appendage be attached only to a rigid body, or to more 
than one rigid body if those bodies are themselves inter- 
connected in such a way that their relative motion cannot 
induce deformation of the flexible appendage. Thus the 
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appendage cannot properly be the only connecting struc- 
ture between two rigid bodies, but it can be hinged to two 
bodies if their relative motion is restricted so as to permit 
die appendage to accommodate their relative motion by 
rotating without deforming. Figure 3 illustrates several 
examples of idealized mechanical assemblages that might 
present real or imagined obstacles in modeling for hybrid- 
coordinate analysis. In each of diese sketches the struc- 
ture drawn as a truss or gridwork is supposed to be sub- 
stantially more flexible than the solid components. 

Figure 3a bears some resemblance to a dual-spin vehicle 
with a flexible antenna, two flexible solar-cell arrays, and 
two flexible booms on die rotor. Here the de-spun plat- 
form Bi, the nutation damper mass B 2 , the rotor B 3 , and 
the contents of a spherical tank of fluid B are idealized as 
rigid. Substructures A 1 , ■ ■ ■ ,A r ‘ are idea ized as elastic 
appendages. An arbitrary control system determines the 
relative rotation of B, and B ;l , with the use of data from 
sensors anywhere on the vehicle. In addition, the entire 
antenna structure A 1 may be controlled to rotate on die 
base B t without exceeding the capacity of the method. 
But if the dish portion of A 1 is controlled to rotate relative 
to its support tower, the model must be modified if the 
equations to follow are to be applied. An individual 
appendage is permitted only elastic deformations, so 
antenna dish and tower could not together be one append- 
age if large relative rotations are prescribed. Yet they 
cannot be treated as two separate flexible appendages, 
since each appendage must be attached to a rigid body 
(and not to another appendage). The only recourse is to 
model the tower as one or more rigid bodies, in which 
case the derivations of the following sections continue 
to apply. 

Similarly, in Fig. 3b one must treat the two flexible 
antenna structures, together with their interconnecting 
rigid body, as one flexible appendage (as indicated by 
the dashed lines). A given flexible appendage cannot be 
attached to two rigid bodies for which relative motion 
deforms the appendage, so the massive block between 
the two antenna dishes cannot be isolated as a rigid body, 
and must instead be absorbed as part of the larger 
appendage A. 

Again in Fig. 3c one is strictly precluded from treating 
each of the three blocks in the system as a rigid body; the 
block shown within the dashed lines must be included as 
a part of the appendage. Because of the large relative 
motions permitted between B 1 and B 2 , one does not have 
the option of including B x plus B 2 with the flexible struc- 
ture instead of the block indicated. If an articulated ele- 


ment such as Bo were attached as well to the block shown 
within the dashed lines, the hybrid-coordinate method as 
developed in this report would be applicable only by 
ignoring the masses on one end of the appendage and, 
instead, applying appropriate forces and torques to the 
end. This represents some distortion of the method, 
and jeopardizes the rationale to be applied to modal- 
coordinate truncation. 

Figure 3d illustrates one of the ways in which a flexible 
appendage can be attached to two rigid bodies. Con- 
straints between bodies B x and B 2 are such that their 
relative motion results only in the rotation (not in the 
deformation) of the appendage A. It would not be desir- 
able in this case simply to include the little body B-. as 
part of the flexible appendage because of the discrete 
damper connecting B 2 to B x . The attachment between a 
flexible appendage and a rigid body must be consistent 
with the assumption of modal damping of the appendage 
if modal coordinates are to be obtained for second-order 
equations. 

When two different points of a flexible structure are 
attached by discrete dampers (e.g., dashpots) to a rather 
rigid base (as in Fig. 3e), one might assume that an addi- 
tional rigid body (such as B 2 in Fig. 3e) connects the 
support points. If this assumption is unacceptable, and a 
discrete damping mechanism must be incorporated within 
the appendage, then one can find modal coordinates only 
after rewriting the equations of motion as first-order equa- 
tions. Although modal transformations of first-order equa- 
tions with arbitrary damping are treated in Section III-D, 
it is tentatively assumed in Section III-B that any damp- 
ing in the appendage may be modeled as classic modal 
damping (Ref. 30). 

It is evident that the equations in this section will not 
apply to completely general models of space vehicles, and 
even when applicable they may not constitute the opti- 
mum approach to simulation. For the system of Fig. 3c, 
for example, it may be preferable to ignore entirely the 
mass in the truss structure and to treat the vehicle as 
three discrete rigid bodies, perhaps even in a point- 
connected topological tree configuration if the truss is 
short and longitudinally stiff. A discrete-parameter simu- 
lation of the vehicle in Fig. 3a would, on the other hand, 
be hopelessly inefficient, and the results would be inferior 
to a much less expensive hybrid-coordinate simulation 

The traditional practice of normal-mode-coordinate 
transformation of the second-order equations of motion of 
the entire vehicle could not be applied directly to any of 


14 


JPL TECHNICAL REPORT 32-7329 



RIGID BODY B 
SOLAR PANEL A 3 













the vehicles of Fig. 3, although for restricted modes of 
motion one could obtain vehicle normal coordinates for 
the corresponding first-order state equations. This pro- 
cedure might, for example, prove optimum for restricted 
motions of the mechanisms shown in Figs. 3d and 3e. 

In reading the present section it may be helpful to keep 
Fig. 3a in mind, since this is the kind of system for which 
the hybrid-coordinate approach is most ideally suited. 

In what follows, attention is focused first upon an indi- 
vidual flexible appendage A attached to a single base B. 
In the exceptional case (e.g., Fig. 3d), where a second 
rigid body is attached to the appendage, the influence 
of this body is felt in changing the attitude of the pri- 
mary base and the appendage in its undeformed state. 
This attitude is reflected in a direction cosine matrix C 
which is considered time-dependent in the derivatio 

After the individual appendage equations are exam- 
ined in detail, the equations of motion of the balance 
of the vehicle are considered. There is no attempt to 
include in this report any general procedure for approach- 
ing the derivation of these equations. A class of vehicle 
somewhat less general than that illustrated in Fig. 3a is 
treated here in detail, although for more complex systems 
there may be some advantage in combining some of the 
discrete-parameter methods of Section II with the ap- 
pendage equations derived explicitly here. 


mechanisms are excluded from th*s idealization,'* but 
modal damping is incorporated in the equations at a 
later stage of the derivation, with the introduction of 
modal coordinates for the appendage. At this point, the 
equations will lose their apparent restriction to discretely 
modeled appendages, and will be equally applicable to 
continuously modeled elastic appendages. In special ap- 
plications to taut strings and membranes or to uniform 
beams, plates, or shells, the analyst may find a continu- 
ous model most convenient. In the vast majority of space 
applications, however, a discrete-parameter model is 
necessary. Frequently the model consists simply of elas- 
tically interconnected particles, but the present assump- 
tion of interconnected rigid bodies is more general and 
often more convenient. 

Let the set of dextral orthogonal unit vectors b,,b^, b ;i 
be fixed in the base B, and the similar set a,, a.., a ., be 
fixed in . e reference established by A prior to defor- 
mation (i.e., fixed in that portion of A contiguous to rigid 
body B). When A is elastically connected to B, these sets 
of vectors may be identical, and at most are related by 
a constant transformation matrix. But when A can rotate 
relative to B (as in Fig. 3a, the antenna A 1 might rotate 
relative to the base B,), the transformation matrix relat- 
ing these vector bases will vary with time according to 
an independently specified control law. If the unit vectors 
are written as vector arrays 


(a t 

1 

( 

{a}- a. 

•, (M — j 

/ 

1 a. 

( 


B. Flexible-Appendage Equations 

Consider at the outset the equations of motion of a 
linearly elastic structure A attached to a base B, assum- 
ing that the structure undergoes only “small” defor- 
mations while the base motion is arbitrary. (“Small'’ 
deformations means in a strict mathematical sense arbi- 
trarily small or vanishing deformations, since terms above 
the first degree in scalar measures of deformation and 
deformation time derivatives are ignored completely. In 
engineering practice, however, this requirement is inter- 
preted quite liberally, and a beam that suffers a tip de- 
flection of 10$ of its length is commonly described as 
having “small” deformations.) 

For convenience in derivation and for compatibility 
with engineering practice, the flexible appendage A is 
idealized initially as a collection of elastically intercon- 
nected, discrete rigid subbodies A u • • * ,A„. Damping 


these arrays may be treated like column matrices in re- 
cording the relationship between {a} and {b} in items of 
the direction cosine matrix** C, i.e., 

{a} ~- C {b} (34) 

In application it is often convenient to be able to select 
the basis (a) for each appendage individually, and to 
select the basis {b} independently (perhaps guided by 
the principal axes of inertia of the total vehicle, or by 


‘‘•Terms to Ik* added to the equations of this .subsection to accom- 
modate idealized discrete dampers are easily generated, and are 
accommodated in the coordinate transformations of subsequent 
subsections. 


•"“in this report matrices are not identified by brackets [ 3 unless 
the elements of the matrix are recorded in detail. Vector arrays 
are always enclosed in braces { }. 
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design coordinates of a drawing, or by attitude-control 
axes). Thus even when A does not rotate relative to B, 
the matrix C may not be the identity matrix E. When A 
does rotate relative to B, so that C varies with time, it 
will become necessary either to augment the equations 
in this section by added dynamic equations or to assume 
that perfect control results in known time behavior of C. 

In devising the discrete-parameter model of A, it is 
convenient to select basis {a} and the individual mass- 
center principal axis bases of A,, • • • , A h to be identical. 

The first step is the formulation of the equations of 
motion of a typical body A, s . in form convenient for com- 
bination with corresponding equations for all n bodies. 
Consequently the vector-dyadic equations of translation 
and rotation will be written as matrix equations in 
basis {a}. 

Let P* be the mass center of body A.,, and let O' be 
a point fixed in inertial space (see Fig. 4). Then Newton’s 
second law provides 


F* - m, (P,) (35) 

where F 4 is the resultant force applied to A*, tn x is the 
mass of A s , P s is the vector from O' to P x , and (as noted 
previously) each dot over a vector denotes time differ- 
entiation in an inertial frame of reference. 

As shown in Fig. 4, the inertial position vector P,, may 
be written as 


P* = X + c + R + r* + u s (36) 


Since a matrix formulation is ultimately required, 3 by 1 
matrices are defined for each of the vectors in Eq. (36) 
in terms of the appropriate vector basis. These matrices 
are introduced with the use of vector arrays (as in Eq. 33), 
by the definitions 


X = (i.Li,) 


x/ 

x, 

X, 


{i) r X 


c = (b, bj b :t } 


Ci 

C-. 

C ;t 


“<b } y c 


R — • (bi bj b ;j } 


■R,' 

R, 

R, 


{b} T R 


(37) 


r s ' = {a, a« a.J 


r i 

rf, 

f 


3 J 


(a } T r s 


u s = (a t a. a : ,} 


u* 




--■2 {a ) T u * 


where the superscript T denotes the transpose of a matrix 
or vector array. 


where X is the inertial position vector of the vehicle 
mass center CM; e is the vector from CM to the point O 
fixed relative to body B and coincident with CM when 
the vehicle is in some nominal undeformed configura- 
tion; R is the vector from O to an arbitrary point Q fixed 
in B on the interface between A and B; r* is the vector 
locating from Q the point (),, occupied by P x when A 
is undeformed; and u* defines the translational deforma- 
tion of the appendage at point Q s . The vectors in Eq. (36) 
are not all expressed conveniently in any one vector basis. 
Vector X defines the vehicle trajectory, which may be 
known in terms of an inertial reference. It thus becomes 
desirable to express X in terms of an inertially fixed vec- 
tor basis ii,ij,i;i. Vectors R and r 4 , on the other hand, 
are fixed in the reference frames in which bj,b 2 , b :( and 
ai.aj.a;, are embedded, respectively, so they are most 
efficiently expressed in their own natural vector bases. 


Upon substituting Eq. (36) into Eq. (35), one confronts 
the necessity of differentiating in an inertial reference 
frame a number of vectors that are much more easily 
expressed (and sometimes constant) in some other ref- 
erence frame. This is most easily accomplished by liberal 
use of the identity from vector differential calculus 

hr 1 hr] 

^ V= _ V + ^ X V (38) 

where V is any vector, is the angular velocity of 
any reference frame relative to any other reference 
frame f u and the superscript preceding the derivative 
operator denotes the reference frame of differentiation. 
Although t may be any independent variable, in the 
present context it is always time. 
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Fig. 4. Discrete-parameter appendage 
sub-body coordinates 

If the presuperscript i denotes the inertial reference 
frame, Eqs. (35) and (36) may be combined as 


where b is the reference frame of body B and basis {b}, 
and co is the inertial angular velocity of b. Differentiating 
Eq. (41) provides (with the use of Eq. 38 again} 

i c p b ( jt 

IF (c + R) = IF <° + R) 

b d 

+ 2co X (c + R) + co X [co X (c + R)] 

+ co X (c + R) (42) 

Now Eq. (37) may be utilized, noting R to be constant, 
to obtain 

~ (c + R) = {b} T c + 2co X (b} r c 

+ coX [co X {b} T (c + R)] 

+ cb X {b} r (c 4- R) (43) 

To express this result in matrix terms, one must obtain 
the matrix equivalent of a vector cross product. 

For two arbitrary vectors V and W, expressed in terms 
of an arbitrary vector basis (e), the vector cross product 
may be written as 

VXW= {e} 7 V X {e} r W = {c) r VW 


F s = m. 


W 

clt 


<cl- 


d- 


7 x + dr- ( c + R ) + ^ rS + u *) 


dt- 


where 


(39) 


The first term in brackets in Eq. (39) is (with the use 
of Eq. 37) simply 



" v, l 



V — 

v s 

R 7 - - 

w. 


_v,_ 




gj; ({!}'' X) = (iJ’-X (40) 

since the inertial time derivative of {i} is zero. (Recall 
that a matrix derivative such as X is just the matrix of 
differentiated elements.) 

The second term in brackets in Eq. (39) may be obtained 
(with the use of Eq. 3S) from 


Vs 


o ~v 3 w 

V a 0 -V t 


-V a V, 0 


(44) 


as may be confirmed by expansion. In equivalent terms, 
the matrix product V\V is isomorphic to the vector product 
V X W. Thus the matrix identity 


rv> -u 

VIE = -wv 


(45) 


' d b d 

__ (c + R) ^ _ (c + R) + <0 X (c + R) 


(41) 


follows from the vector identity V X W = — W X V. The 
tilde operator (~) over a 3 by 1 matrix represents the 
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corresponding skew-symmetric 3 by 3 matrix, as in 
Eq. (44), throughout this report. Note the useful identity 

V T =-V (46) 

for any tilded matrix, because of skew symmetry. 


'cl~ a cl 2 °cl 

-j— (r s 4- u s ) = — t— u s + 2 £l a X—r;u s 
dt- ' dt- dt 

+ £l lL X X (r s + u s )] 

i>d 

+ -r-£l a X (r s 4- u s ) 


Equation (43) can now be written (with the definition 
co== {b} r w) as 

(c 4- R) = {b} r [c 4- 2ayc + S£(c + R) + d) (c + R)] 

(47) 

which provides the second term in brackets in Eq. (39). 

The final term in Eq. (39) can be obtained in similar 
form, although the vectors r s and u s are written in basis 
{a} in Eq. (37). The differentiation then is performed in 
two stages. First, 


+ 2wX 


n d 

u s 4- X (r s 4- u s ) 


-}- &) X [to X (r s 4- u s )] 4- tb X (r s 4~ u s ) 

(48) 


where £1" is the angular velocity of the {a} frame rela- 
tive to the {b} frame, and it has been recognized that 
the time derivative of r s in reference a is zero. This result 
is now to be written in terms of the vector array {a} and 
the matrices defined in Eq. (37). If £2" is written in the 
{a} basis as 

Sl a = (a} r ft' 1 (49) 


i d.- h d- b d 

IF (r< + uS) = dr (r * + u-,) + 2co x Ft (rS + l,s) 

+ co X [to X (r* 4- ii s )] 4- co X (r s + u s ) 


and Eq. (34) is used to write 

{b} T = {a} r C 


(50) 


is obtained from Eq. (42), and then the derivatives in so that w is written in the {a} basis as 
{b} are replaced by derivatives in {a} by using Eqs. (41) 

and (42), and replacing i by b and b by a. The result is to = {b} r <o = {a} r Co> 


then Eq. (48) becomes 


1 (J- ^ ^ (N/ 

clF u *) = { a ) Y [w* + 2n" if 4- ft" ft" ( r s 4- if) + ft" (■;•* 4- if) 

4- 2 (Si) if + 2 (Ci) ft" (r s 4- if) + (Cl) (C») (r* 4- if) + (Cl) (r s 4- if)] 




(51) 


When Eqs. (51), (47), and (40) are combined in Eq. (39), the three different vector bases of these equations arc aban- 
doned in favor of basis {a}, noting that this is also the natural basis for F\ This requires the use of Eq. (50) in Eq. (47), 
and the parallel use in Eq, (40) of the relationship 


{i}' r = {b} 7 ’0 = {a} r C0 


(52) 


where (-) is the direction cosine matrix relating the body B to an inertiallv fixed vector basis, i.e., 

{b} = 0 {i} (53) 

The result of these substitutions into Eq. (39) is the vector equation 


{a} 7 F s — {a} 7 ’wj.* {C(r)X H- C [c + 2 tic 4' to(o(c 4* R) 4- <b(c 4- R)] 4- if 4* 2ft" ti s 4~ ft"ft"(r* 4~ if) 

4- ft" (r s + if) 4- 2 (C^) if 4- 2 (6£) ft" (r s - if) 4- (6U) (Go) (r* + if) + (CA) (r» 4- «•)} (54) 
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which employs the new matrix symbol F s as defined by 


F s == {a} 7, F s (55) 

One may now record the corresponding matrix equation (in vector basis {a}), either by inspection of Eq. (54), or by 
formal clot-multiplication by the vector array {a}, noting that the relationship 


la, 1 

1 


'I 

0 

o' 

{a} • {a} r sh |a,j 

► • {a, a 2 i 


0 

1 

0 

/ \ 

\ a <; 

1 


_0 

0 

1_ 


== E 


(56) 


follows from the introduction of the dot between vector arrays to imply the multiplication of vector arrays follow- 
ing the rules of matrix multiplication, but with dot-multiplication of the vector elements inside the arrays. This 
notation is sometimes convenient for notational bookkeeping, but even without this artifice one can simply record the 
matrix equation isomorphic to the vector equation (54). In the process, the future direction of the derivation may be 
anticipated by the replacement of the scalar m s with the matrix m* — m s E, with E the 3 by 3 identity matrix. The 
desired matrix equation is then obtained (after some factoring and rearrangement) in the form 

F* = m s { C [aX + c + 2 Sc + 7SS (c + R)-(c + R) »] - (P -f t*«) (6 fl + Cd>) 

+ [(ft) (&Z) 4- 2 (&) (n«) + n« £2«] (r s + II s ) + 2 [9." + (C^)] w* + «•} (57) 


Of course the matrix c, which denotes the translation 
in body B of the vehicle mass center, depends in part on 
the translations u* (s = 1, • • • , n) of the sub-bodies of 
the appendage A. In general, it may also depend on the 
deformations of other elastic appendages of the system, 
and on the (possibly large) relative motions of other 
articulated rigid bodies of the vehicle. To emphasize the 
influence on c of the deformation of the appendage for 
which equations of motion are being written, the mass- 
center definition is employed to write for c 

1 »' 

c = — -r— 2 m„u # + e (58) 

V'l Sil 

where e accommodates the mass-center motion due to 
other appendages and moving parts, if any exist, and Q ]] 
is the total vehicle mass. With the definition 

cs {b) J ’e (59) 

Eq. (58) may be written in the form 

{b) r c = — — yy 2 {a } T m*u* + {b) T e 

= {b } r [-2CV«* + e] 

8 - 1 

where Qll and Eq. (34) has been employed. This 

vector equation has the matrix counterpart 


c — — 2 CVu* + e (60) 

S = 1 

Thus the role of the appendage deformation in Eq. (57) 
is made more explicit by the substitution of Eq. (60) and 
its derivatives (permitting C to vary with time), which 
results in 

6 - - 2 (CVu* + CViP ) + e (61) 

S - 1 

and 

c = -2 (CV«* + 2 CVu* + C r iAi s ) + e (62) 

If it seems preferable to avoid the time derivatives of the 
direction cosine matrix C, this can be accomplished by 
using the identity (38) once more to obtain 

{b} r c = ^[{b}'c3 

= 4; [-w 2 (* s «' + {b}’ e] 

* S ; ; 1 

+ X [ — {a} 1 ' 2 /**«*] + {b} 7 ’ 6 
1 

= —{a} r 2 — {a) r Q ff 2 + {b} r e 

Sl s-l 

(63) 

t • ✓ 

JPL TECHNICAL REPORT 32-1329 


20 



With Eq. (34), this yields the matrix equation 


c = —2 CV« S - c ’ r 3 a 2 i-^u s + e 


Differentiation of Eq. (63) in reference frame b similarly provides the vector equation 


(64) 


(b} r c = — {a} 5 " [ 2 /Wi 8 + 2 n 3 u s + Q a 2 P 3 u a ] - (a) 7 fi" [ 2 n a u s + 0“ 2 /**»•] + {b} r e 


8 = 1 


S = X 


The matrix counterpart in basis (b) for this vector equation is 


c = -C T [ 2 M s « 3 + n a 2 + 2n a 2 p a u a + 3 ,l n a 2 |U 8 W S ] + e 


(65) 


Now Eqs. (64) and (65) may be used in Eq. (57), instead 
of substituting Eqs. (61) and (62) in this equation, thereby 
avoiding derivatives of C. Alternatively, one can avoid the 
relative angular rate Q' 1 and its time derivatives entirely 
by the substitution of the kinematic identity (Ref. 16) 

n fi = CC T (66) 

and its consequence 

n* = CC T + C C T (67) 


With these substitutions in Eq. (57), and with Eqs. (61) 
and (62), all relative motion between bases {a} and (b) 
is expressed in terms of the direction cosine matrix C 
and its derivatives. Although this option may be optimum 
in some applications, it would seem preferable in most 
practical cases to avoid C and C, and therefore to employ 
Eqs. (57), (64), and (65). 


Since the sub-bodies A x . • • • , A„ are considered to be 
rigid bodies rather than particles, the translational equa- 
tion (57) must be augmented by rotational equations 
such as 

T s = H 8 (68) 


for typical rigid sub-body A s . Here H s is the inertial an- 
gular momentum of A,, referred to the mass center P s , 
and T s is the corresponding resultant torque. The rules 
of vector differentiation (see Eq. 38) may now be applied 
to H 8 , or to its equivalent I s • co 8 , where I* is the mass- 
center inertia dyadic of A s and co 8 the inertial angular 
velocity of A s . The result is 


— H s ) 
dt 

-r (i* • oA 1* • — & > 8 + -f i 8 • to* 

dt ' dt dt 


hi 

= I 8 • to 8 -f- 
dt 



I 8 + co 8 X I 8 


- 5 s X <o 8 ^ • co 8 

(69) 


where the symbol s d/dt denotes differentiation in the 
reference frame of A s . Here use has been made of the 
counterpart for dyadic differentiation of the vector dif- 
ferentiation relationship of Eq. (38), namely 


<hI 

dt 


D = 


hd 

T D + X D - D X '>co'» 
dt 


(70) 


with D any dyadic and f x and f» any two frames of ref- 
erence. (This relationship can be confirmed by writing 
the dyadic in the form* D = D a p e« ep, and then applying 
Eq. (38) to the unit vectors e« and ep.) 

The term I s X co 8 • to 8 in Eq. (69) is zero, since the 
implied operations include dot-multiplying by co 8 a vector 
orthogonal to co 8 . Equation (69) is further simplified by 
the constancy of the inertia dyadic of the rigid body A s 
in the reference frame of that body, i.e., ( *d/dt ) I 8 is 
zero. In rewriting the correspondingly simplified rota- 
tional equation, one may expand the inertial angular 
velocity co 8 of A s by use of the “chain rule” 




(71) 


"Throughout this report, lower-case Greek indices range in value 
from 1 to 3, and when these indices are repeated in a given term, 
summation over these values is implied. 
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where • • • ,f n are n arbitrary frames of reference and 
f ’d fr is the angular velocity of frame f r relative to frame 
In this case the chain is from the frame of body A* to 
the frame of A, to the frame of B, to inertial space. Since 
the rotation of A., relative to A is due only to “small” 
structural deformations, this rotation can be represented 
by the vector p s = /3 8 a t + /3? a 2 + /3 8 a 3 , where /3J./3!!,/? 8 
are three angles of rotation about oithogonal axes ai,a 2 , a 3 . 
This representation is strictly applicable only if these 
rotations are infinitesimal, but this is nonetheless the com- 
monplace assumption of structural dynamicists. Conse- 
quently co s in Eq. (69) can be replaced by the use of 
Eq. (71) to write 

£o 8 = CO + S2." + {a} 7 /3* (72) 


where /3 8 is the 3 by 1 matrix defined by fi s = (a } T (3 s . 
This substitution puts Eq. (69) in the form 

'cl 

T 8 = I* • (CD + SI" + {a} r (3 s ) 

+ (<a + SI" + {a} 7 /3 s ) X I s • (to + S2." + {a} 7 /3 s ) 

(73) 

Again the differentiation is accomplished with the help 
of Eq. (38), and simplification is afforded by the assumed 
smallness of /3 s (since only first-degree terms in /3 s and 
f3 s are retained). The result may be written as 


•p — 


'cl 

clt 


03 


"d 

dt 


SI" -f (c o + SI") X SI" + {a} 7 ’ /3 s + (« + SI") 



•f (a + SI") X I 8 • (w + SI") +• (&> + SI") X I* • {a} 7 ’ /3 s + {a} 7 j3° X I* • (« + SI") 


(74) 


If all vectors are now written in terms of the vector arrays of the most convenient basis, with the use of the defini- 
tions following Eq. (49) together with 


T 8 s={a} 7 T 


and, with the summation convention, 


I 8 == I s afs a n af} == {a,a s a s } 


— 1 

I s 

1 12 

I S 

I 8 

L 31 

I s 

■*32 


I s ' 

x ia 


I s 

A 113 

I S 

*33 , 


at 

a 2 ) {a} 7 I s {a} 

a* 


the vector rotational equation becomes 


{a} 7, T a = {a} T I s {a} • [{b} 7 'J> + {a} 7 6“ + {b} 7 * X {a} 7 Q“ + {a)' r (3« + ({b} 7 w + {a} 7 n") X {a} 7 /3 s ] 

+ ({b} 7 <o + (a} 7 fi a ) X {a} 7 / 8 {a} • ({b} 7 «. + {a} 7 SI") 

+ ({b} 7 « + (a} 7 n fl ) X {a} 7 I s {a} -{a} 7 /3 8 + {a} 7 /3 8 X {a} 7 ! 8 {a} *({b} 7 o> + {a} 7 0°) 


(75) 


(76) 


By using Eq. (50) to obtain every vector in basis (a), and applying the matrix representations of vector cross- 
multiplication (see Eq. 44) and vector array dot-multiplication (Eq. 56), one can obtain this vector equation in the 
useful form 

(a} 7 T 8 = {a} 7 {7 8 [CA + SI" + (62) Sl a + f3° + (62) /3 8 + Q a /3 8 ] 

+ [(62) + fi n ] I s [Co, + X) a ] + [(62) + n a ] P/3 S + ji 8 / 8 [Cm + Q a ]} (77) 

Now the isomorphic matrix equation can be written by inspection (or obtained formally by dot-multiplying by {a}). 
In the process, the identity of Eq. (45) is applied whenever it seems convenient to have the most obviously unknown 
variable available for factoring on the extreme right. For example, in the last term the unknown deformation /3 s must 
appear on the right if it is to be factored out with the term preceding it, so the identity 

/3 8 I s (Ca> + SI") = - [I s (C<0 + fi a )]~/3 8 (78) 


22 


JPL TECHNICAL REPORT 32-1329 



is used. Here the tilde outside the square brackets is to be construed to apply to the entire 3 by 1 column matrix within 
those brackets. With this substitution, the matrix equation corresponding to Eq. (77) becomes* 

T' = p [Ci + n» + /3 s ] + [I s (C u ) + (Co,) I s + I s o fl + n n I s - ( I s C M + I s n«) ] /? 5 

- I s n u C<o + [(C^) + S“] P [(Cm) + p. a ] (79) 

Again it should be noted that these rotational equations can, with the use of Eqs. (66) and (67), be written in a form 
not involving fi a , should this be desirable in a special case. 


In most applications, it is convenient to accept Eqs. (79) and (57) as the rotational and translational equations of 
motion, respectively, of body A„, with Eqs. (60), (64), and (65) substituted into Eq. (57) to accommodate vehicle mass- 
center motion relative to body B. The result of this substitution is recorded as follows, with the temporary convention 
that 25 means the sum over k ranging from 1 to n: 


F s — m 3 {CeX — (CR + r s C + Ce) m — r 3 fi° + i/ s — i l k + C'e + 2 [Q“ + (C<o)] ii 3 — 2CS C T 2ji k u k *t* 2 Cm e 

- u k - C%C T 2p*tfi + [o n + (Ci)] u s - n a Xn s u s - 2CsC T n a - n a n Q u h 

+ C SS (R + e) — C SSC 7 ’ u k + [(Cm) (Cm) + 2 (Cm) (Q a ) + (r 3 + u s )} (80) 




Again the identity (Cm) s=CmC T might be substituted, and also the consequence (Cm) (Cm)==CmmC T . 


Equations (79) and (80) constitute the most general 
formulation of equations of motion of an appendage sub- 
body A„ to be derived in this report. These equations are 
applicable to a sub-body of an appendage that is under- 
going small deformations while rotating in an arbitrary 
way (as described by C and Q a ) relative to a base B sul> 
ject to any translation and rotation (as described by X 
for the vehicle CM and 0 and m for body B). Any number 
of additional appendages or articulated moving parts may 
also be attached to B; these internal motions will influence 
0, o), and e. 

In practice, the generality of these equations is rarely 
required. Only in the exceptional case (e.g., the scanning 
antenna) is there a nonzero fi", so C is usually constant. 
Then there is no theoretical objection to replacing C by 
the identity matrix E, particularly if there is but one 
appendage attached to body B. For reasons not yet ap- 
parent, it may be computationally desirable to adopt the 
view that all bodies B have only one appendage, even 
though the appendage may in some cases be composed 
of several physically distinct structures attached to B. 
Then e is zero, except when a statically unbalanced rotat- 
ing rigid body is also attached to B, or a secondary rigid 
body can translate relative to B. 

For the simplest configuration, with C = E and e==0, 
the equations of motion (79) and (80) adopt the form 

T 1 = I s [m + (3 s ] + [I S S + ml* - (Pi) ] /?* + SPm 

(81) 


F s = m s [0X + ii s — 2jx k u k + 2m(ii s — %f.i k ii k ) 

— (R + r s )S + < I) (u s — 2fi k u k ) 

+ mm (R + r 3 + w s - V u k ) ] (82) 

These equations still permit unrestricted motion of base B. 

The forces F s and torques T s that appear in the basic 
differential equations (79) and (80) of sub-body A s include 
terms due to structural interactions with neighboring sub- 
bodies of the appendage A. These interactions are repre- 
sented here as linearly elastic and viscous forces and 
torques, so they are proportional to the deformations 
and deformation rates as represented by u\ * • • , u", 
f3\ • • • , /?" and their first derivatives. In practice, viscous 
damping terms are often unspecified until transformation 
to modal coordinates has been accomplished. 

The resulting approximation is far from perfect in its 
representation of structural behavior. Even when the 
materials of a complex structure are essentially linearly 
elastic, there may be sufficient “play” in the joints of the 
assembled structure to jeopardize the assumption of linear 
elasticity of the composite structure. Furthermore, the 
assumption of viscous damping is not strictly in conform- 
ance with the performance of even the simplest structural 
elements. But for complex structures lacking discrete arti- 
ficial dampers, the assumptions of linear elasticity and 

*The identity (6«) = CmC 1 ' should perhaps be noted, although 
this substitution into Eq. (79) seems computationally inefficient. 
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viscous modal clamping are almost universal practice in 
structural dynamics, and these assumptions are reason- 
ably well supported by vibration test data (Refs. 31, 32). 
They are accordingly the assumptions adopted in this 
report. 

Structural forces and torques in each of the equations 
typified by Eqs. (79) and (SO) therefore may couple these 
equations to those of every other sub-body of the ap- 
pendage. The most convenient method of recording all 
the coupled translational and rotational equations of the 
n bodies A,, • • • ,A„ is with a single matrix equation of 
dimension 6n. 

Because interest is focused presently on the appendage 
deformations « l , • • - , ?/", /3\ • • • ,/?“, only terms in 
these variables will be written on the left side of the 
appendage matrix equation. The single 6 n by 1 matrix 

q^[u\ til til P) PI PI ll 'i ■ ■ ■ P'iY (S3) 

fully characterizes the appendage deformations. 

Because the equations have been linearized in the ap- 
pendage deformations, the matrix equation must have the 
structure 


nition, since the matrix dimension is always clear from 
the context in which it is used. 

The operators Seo and S„e serve to distribute a given 
3-row matrix into alternate 3-row partitions of a larger 
matrix, having perhaps 6n rows. For example, the expres- 
sion C(-)X is the column matrix 

[(C0X) 0 (C(-)X) 0 • • • (C0X) 0] 7 ' 

which, by inspection of Eqs. (79) and (SO), must appear 
as part of 1/ in Eq. (84). 

The operators % T E() and are, on the other hand, 
effectively selective summing operators, since, for exam- 
ple, when multiplied by a 6n-row matrix they sum respec- 
tively the odd- and even-numbered partitioned 3-row 
submatrices of the 6m -row matrix into a single 3-row 
matrix. As a specific example, q provides the 3 by 1 
column matrix obtained by adding u 1 + ir + • • • + ?/". 

The constant matrix M f in Eq. (84) is now most easily 
written in terms of the new symbol M, defined in terms 
of 3 by 3 partitioned matrices as 


A/' q + D' q + G' q + K' q + A' q — U (S4) 

where by definition D' and K' are symmetric and G' and 
A' are skew-symmetric. Inspection of Eqs. (79) and (SO) 
reveals that A/' is a constant symmetric matrix, but G', 
A', and K' depend on the variables <o, Q", and C. The 
matrix U depends on these variables and in addition on 
0, X, e, and the external applied forces that may appear 
in F* and T*. The matrix D' accommodates damping in 
the structure. 


m 1 0 0 0 • • • 0 

0 7 1 0 0 • • • 0 

0 0 m- 0 • • • 0 

M = 

0 0 0 P • • • 0 


0 0 0 0 • • • I* 


Now M' can be recorded as 


( 86 ) 


The detailed representation of the 6m. by 6m coefficient 
matrices A/', G', A', and K' and the 6m by 1 matrix U is 
facilitated by the introduction of the Boolean operator 
matrices 

S so = [E 0 E 0 • ■ • E OF ) 

> (85) 

= [0 E 0 E ■ ■ ■ 0 E] T ) 

where E and 0 are 3 by 3 matrices, the former being the 
unit matrix and the latter the null matrix. In the present 
application the matrices Seo and % 0E will generally have 
6m rows, but this restriction is not embodied in the defi- 


A/' = M (E - 2 £0 2L M/Qli) (87) 

Here E is an identity matrix of dimension 6m by 6m, 
and 111 is, as previously, the total vehicle mass (recall 
ix s s= m >l /7}li). The matrix M may be recognized as the 
inertia matrix that would be involved in a calculation of 
cantilever modes of the appendage. This interpretation 
will be explored in detail in Section III-D on coordinate 
transformations. 

In expansion, A/' appears in terms of 3 by 3 partitioned 
matrices as shown in Fig. 5. Note that M' is a symmetric 
matrix of constants. 
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m x — mW/Q?? 0 —mtmr/Qft 0 


0 


0 0 


0 m- — m z m 2 /(Y)l 0 


0 0 0 7- 


M' = 


—m n m 1 /Q')i 0 —m n m 2 /Q7i 0 


0 0 0 0 


Fig. 5. Inertia matrix 
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Before recording the matrix G' in expanded form, it is useful to examine the coefficient of {3* in Eq. (79) in some 
detail. It is evident that neither Q" nor C changes the character or structure of this coefficient, so the examination is 
restricted to the special case for which C = E and Q" = 0. By expanding the matrices, one finds that 


0 — (E 4- Z 2 — Z 3 ) q) 3 

Z*S> + SI* - (K) = (7 X -r L — l.) o J: . 0 

(7: - Jo + I a ) co 2 ( — Zj -r L - 7,) Wl 

= (tr 7 s ) 5 — 2 (TV) 


Thus it is evident that the sum of these three matrices 
is skew-symmetric, although tire first two matrices are not. 
With Eq. (SS) substituted into Eq. (79), inspection of this 
equation and Eq. (SO) permits the matrix G' to be recorded 
in terms of 3 by 3 partitioned matrices in the expanded 
form shown in Fig. 6. The substitution C<SC T = (Co) has 
also been applied here. Note that all of the q terms on 
‘the right side of Eqs. (79) and (SO) are accommodated 
in G', leaving only damping terms in the symmetric 
matrix D'. 

The matrix G' is skew-symmetric, since the transpose 
of the matrix in Fig. 6 is also its negative. A more compact 
and useful representation of this matrix can be obtained 
by broadening the definition of the operator tilde (~). As 
illustrated in Eq. (44), when a tilde appears over a symbol 


signifying a 3 by 1 matrix, the corresponding skew- 
symmetric 3 by 3 matrix is implied. Henceforth, when a 
tilde appears over any symbol representing a column 
matrix that can be partitioned into 3 by 1 matrices, a 
skew-symmetric matrix is implied, which, as a matrix of 
3 by 3 partitions, is diagonal, with the 3 by 3 skew- 
symmetric counterparts of the 3 by 1 elements of the col- 
umn matrix ranged along the diagonal. For example, on 
the column matrix 


O' 1 = [Q“ 0 Q° 0 • * • Q* 0] 


the tilde oyoerator signifies 
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With this added convention, after inspection of Fig. 6 and Eqs. (79) and (SO), the matrix G' may be written in the 
compact form 


G' - 2 M {(S f „flT+ (2* tl Co,)~- S* u [(Co,) + Q"] M/072]} + M [(2„*a«)~-l- (2o*Co,)~] 

+ [( 2 „*Q“)~ + (s 0 J! C«)~] M - [M (% 0E a a + ZobCco)]' 


(90) 


The matrices A' and K' appear in compact form as 


A' = M {[(2*„n*r+ (s*„Cin - 2™ [O' 1 + (Ci)] tU M/Q71]} \ 

K' = M {(S*„ O")~(2* 0 £2°)~- + 2 (2*o Ca>)~(2*„ 0“)~+ (S«o Co>r(S*„ Co ,)~ f (91) 

- 2*„ [£2» Q n + 2 (&) O' 1 + (Ct) (&)J S* 0 Af/07Z} + K ) 


where K is the stiffness matrix that defines the structural 
interaction forces and torques induced by deformation of 
the appendage. The stiffness matrix K is symmetric (see 
Refs. 11-13), but in general there does exist a skew- 
symmetric matrix A', Specifically, the terms in Eq. (91) 
involving angular accelerations 'V' and d, are skew- 
symmetric, and all other terms are symmetric. 

Figure 7 shows the first few terms in the upper left 
corner of the matrix A' 4 K' — K. It is important to rec- 
ognize that in checking this matrix for symmetry one 
must, in transposing the partitioned matrix, remember 
to transpose the 3 by 3 matrices within the partitions. 
Since the transpose of a skew-symmetric matrix is its 
negative, one has, for example, 

(n a ) T = -a a 

and 

r rv> nJ _ ... ,no . ... .no no . . no . no no 

[n°n B ] r = ( n a )‘ (a 11 ) 1 = (-n a )(-a a ) = ci"n a 


This permits the previously noted identification of the 
angular acceleration terms as the skew-symmetric terms 
that compose A'. 

There remains only the explicit identification of the 
elements in the matrix U of Eq. (84). Compact represen- 
tation is accomplished here with the new definition 

/• === [r 1 0 r- 0 • • • r" 0]'' 1 (92) 

The general interpretation of the tilde operator now 
provides 



0 

0 

0 

... 0 

0" 

0 

0 

0 

0 

... 0 

0 

0 

0 

7 - 

0 

... 0 

0 

0 

0 

0 

0 

... 0 

0 

0 

0 

0 

0 

, , . 

0 

_0 

0 

0 

0 

... 0 

0_ 


Now the 6n by 1 matrix IJ can be written as 

V - - M 2 n * [a a + Cm - Q a Ct 9) “ [2oe + Co,)] ~ M 2o* (Q“ -I Co,) 

- M 2*o C [®X - (R + e) 0 , -!- e - 2£o, + SS (R + e)] 


+ M {r 2* 0 (n« + Cd,) - [(2*oWr(2* Q n'‘)^+ 2(S JJ0 Ca>r(2 fl , I fl«)' - + (S#„ C«f (S*„ C<o)~] r) + A 


(94) 
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where the new symbol A is a column matrix in which to 
store any forces or torques applied to the sub-bodies 
A,, • • • , A„ other than the structural interaction forces 
induced by deformations. For example, field forces such 
as gravity forces on the appendage belong in A, and any 
attitude-control jets on the appendage would contribute 
to A in a manner established by a control law. 

In expanded form, the matrix JJ appears as a column 
of 3 by 1 submatrices in Fig. S, in which the forces and 
torques on body j in A are represented respectively by 
f j and V. 

This concludes the derivation of the equations of motion 
of an elastic appendage on a base undergoing arbitrary 


motion. The final result is Eq. (84), with matrices q , 

Q\ K\ A', and L' given in compact form by Eqs. (83), (87), 
(90), (91), and (94). These equations are complicated by 
the generality introduced with time-varying C, which 
permits application to a flexible appendage of changing 
attitude relative to its base (e.g., a scanning antenna). 
Additional complexity stems from the variable e, which 
is present to accommodate movable parts within the ve- 
hicle (other than the appendage). Specialization to the 
simplest configuration, with Cssfi and e 0, led to the 
basic Eqs. (81) and (82) for the typical submass A,,. 

When these equations are combined as a single matrix 
equation such as Eq, (84), one obtains the simpler result 


M (E - Z Bn M/qn) q 4- {2M [(S*, <.,)~ - 2» S 2 *, M /<)]}] 4 - At ( 2 „« 4 (S„* «>)~M - [M v„« « ]~ I D'} q 


+ {M [(2b« w)~ — m , M/Qli] 4- M M/Qll] 4- K) q 


A1 22 1 1 /. * o> (Ai>; v ' to) A'l (Aiuj w) ' A'l ( A/^(j [ (-) A Rw 4* o> ii> H] r 22 /.; ( , <d (2b» w) (2/jn «>) r) 4 A 


(95) 


Even in this case the equations remain complex in 
appearance. Although the specialization has greatly re- 
duced the number of terms in the equation, no simplifica- 
tion of the basic structure of the equations has resulted. 
For the special configuration as well as the more general 
one, the equations of motion have the structure of Eq. (84). 
In both cases Al' is constant and symmetric; G' is skew- 
symmetric and depending linearly on the variable and 
generally unknown angular velocity w of the base 73; 
K' is the sum of a constant symmetric stiffness matrix K 
and a variable symmetric matrix depending on second- 
degree terms in <»; A' is a variable skew-symmetric matrix 
depending linearly on «>; and 7/ is composed of applied 
forces and torques, nonlinear functions of «, linear terms 
in (1), and the unknown mass-center motion X. 

It may appear that the various matrix operators intro- 
duced here serve to obfuscate rather than clarify results, 
and that the appendage equations are too hopelessly large 
and complex for efficient simulation. Both of these impres- 
sions would be correct if the derivation were to end at this 
point. Many of the matrices appearing here are quite 


sparse, and their use in a digital computer program is 
rendered inefficient by the number of multiplications by 
zero. The number of sub-bodies Aj, • • - , A„ required for 
an acceptable simulation of a flexible appendage may be 
quite large (30 is not an unreasonably large value for n), 
and the matrix Eq. (84) is equivalent to 6n simultaneous 
coupled scalar differential equations. 

Practical utility can be demonstrated for these equa- 
tions only after transformation to carefully selected 
coordinates permits the exercise of judgment in ignoring 
all but a small number of variables. It is remarkable that 
in introducing a coordinate transformation and combining 
the resulting equations with appropriate equations for the 
balance of the vehicle, seemingly complex matrix combi- 
nations will emerge with simple physical interpretations, 
and coordinate truncation will permit meaningful prelim- 
inary designs of vehicle attitude-control systems to be 
accomplished without recourse to the computer. The re- 
sults of the coordinate transformations, however, will be 
deferred until Section III-D so that we may first consider 
the remainder of the equations of motion of the vehicle, 
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CO 

o 


A' + K' - K = 


[a a + (62) + n a a a + 2(62) 0 ° + (£2) (62)] [m 1 - mW/)??] 

0 

/y /V ^ ^ ^ /x; rv/ 

— [Q° + (C<i) + fl a Q a 4- 2 (Co,) n a + (Co,) (Co,)] mtmr/Qn 

0 


0 — [Q° + (62) + Q“ O a + 2 (62) + (Cm) (62)] m 2 m\/Qn 

0 0 

0 [6° -f (Co) -f Q“ o a + 2 (Co) Q“ + (Co) (Co)] [to- — m-mr/lli] 

0 0 


Fig. 7. Modified stiffness matrix 


L' = 


■» *• , _ . - . 

— to 1 (C [©X — (R + e) o + e — 2eo + o o (R -re) — r l (Q a + C<i>)] + [ft“ Q a + 2 (Co) + (Co) (Co)] r 1 } -f- / 1 

-I 1 [a a -f Co - a a Co] - [8* + (62)] I 1 [n a + Co] + Z 1 
— m 2 (C [0X — (R + e) o -f e — 2el> + SS (fl + e) - r 2 (f2 a + Co)] + [O a Q a + 2 (62) fl a + (62) (62) r 2 ] } + f 

— I 2 [Q tt + Co — Q a Co] — [o a + (62)] I- [n a + Co] + Z 2 


—I" [6° + Co - n a Co] - [a c + (62)] I" [o a + Co] + Z” 


Fig. 8. Forcing function 



C. Vehicle Equations 

Immediately upon shifting attention from the append- 
age deformations to the time behavior of the total vehicle, 
the analyst confronts a basic question. Should he simply 
isolate the rigid body to which the appendage is attached 
and derive its equations of motion, treating the forces and 
torques applied to it by the appendage as '‘external”? 
Or should he instead write equations of motion for a 
dynamic system that includes both the rigid body and 
the flexible appendage? A parallel question was con- 
sidered in Section II, in which augmented-body methods 
were contrasted with nested-body methods. In the con- 
text of hybrid-coordinate analysis, the answer to this 
question is complicated by the anticipated truncation of 
modal coordinates of the appendage. If the analyst iso- 
lates the body to which the appendage is attached, he 
must be sure that in the truncation process he retains a 
valid representation of the interaction forces between 
base and appendage. Because of this difficulty, it seems 
generally desirable to include the appendage with the 
rigid body in the dynamic system, and this will be 
the method used here. In certain cases, however, when the 
vehicle undergoes large changes in configuration, it may 
be most efficient computationally to treat the appendage 
interaction forces as external forces applied to the rigid 
body. Then some artifice is required (such as the “synthetic 
modes" of Ref. 28) to assure that coordinate truncation 
does not invalidate the interaction force expression, Tin's 
method is explored in Section III-F of this report. 

The immediate task is the derivation of equations of 
motion of a typical vehicle that includes flexible append- 
ages. No attempt will be made here to establish a general 
set of equations applicable to virtually all space vehicles. 
Attention is instead restricted to a class of vehicle that 
may be modeled as a rigid body R, with two flexible 
appendages (A 1 and A"), a balanced rigid symmetric 
rotor 7?.,, and a single rigid body R, attached to R, with 
a single degree of translational freedom. 

The vehicle sketched in Fig, 3a does not quite meet the 
restrictions of the vehicle equations to be derived here; 
the rotor T, of Fig. 3a is not rigid and symmetric, since 
it includes an internal body R, and two flexible append- 
ages A' and A\ Figure 9 better typifies the vehicle to 
which the following equations apply directly, It is an 
idealized tri-spin satellite, designed to maintain spin 
stabilization with the rigid symmetric rotor li M while a 
control system between bodies R i and R, ( maintains the 
rigid platform R, in an earth-pointing attitude, The flex- 
ible antenna A 1 is attached to Ri at a hinge that permits 
limited controlled relative motion, A solar-cell array A' J 


is also attached to R ls with an auxiliary control system 
maintaining the required sun-pointing orientation. A 
linear oscillator B» acts as a nutation damper for the 
satellite; its single degree of translational freedom is in 
a direction transverse to the bearing axis between R, 
and R;., 

Although the vehicle of Fig. 9 is of substantial interest 
in itself— and by specializing to omit components or fix 
joints one can use its equations of motion for a very wide 
range of space vehicles— its derivation is included here 
principally for its illustrative value. The analyst must be 
expected to approach the derivation of equations of mo- 
tion of each space vehicle configuration individually, if 
computational efficiency is to result. The appendage 
equations derived in Section III-B, on the other hand, 
should be applied directly to any flexible appendage 
within the definition offered, and should not require 
rederivation for specific cases. 

One may write in application to any material system 

the basic Newton-Euler equations, F ■ : )//a and T ■ II, 

where 111 is the vehicle mass, F is the resultant applied 
force, and a, H, and T are, respectively, the inertial accel- 
eration, angular momentum, and resultant torque, all 
referred to the vehicle mass center. 

In deriving the rotational equation from 

T II (96) 


TRANSLATIONAL 
DAMPER BODY B, 



EARTH-POINTING ANTENNA A 


HINGE 

RIGID PLATFORM S, 

RIGID SYMMETRIC ROTOR 8, 


SUN-POINTING SOLAR-CELL ARRAY A‘ 

Fig, 9. Tri-spin satellite (cross-sectional view) 
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one simply substitutes the angular momentum definition 


H: 


= J p X p dm 


(97) 


where p is the generic position vector from the vehicle 
mass center to a typical differential mass element of the 
vehicle. 


Thus many of the integrals in Eq. (98) become cither zero 
or a simple vector. The immediate result is 


H — I pXp dm + cX c dll — QUc X (a X c) 
+ j p X (co X p) dm 


Define a new vector p as the generic position vector 
from the point O fixed in body B u and replace p in 
Eq. (91) by p + c (recalling that — c defines the dis- 
placement in B, of the vehicle mass center CM from its 
nominal location at O). The result is then differentiated 
(making use once more of Eq. 38) and the integral re- 
written as follows, employing a convention whereby dif- 
ferentiation in the reference frame of B, is signified by 
an overcircle (e.g., p): 


p X p dm 4- OHc X c T j p X (co X p) dm 

( 102 ) 

The last integral in Eq. (102) may be recognized as the 
dot-product of the inertia dyadic I of the vehicle for point 
O and the inertial angular velocity of the base B,. This 
identification can be established in detail as follows: 


H = 


(p 4- c) X (p 4- c) dm 


— j (p + c) X (p + c) dm 

+ f (p + c) X [co X (p 4- c)] dm 


I p X (co X p) dm — j [p * pco — p • cop] dm 
= j [ p • pE — pp] dm • co I • co (103) 


where E is the unit dyadic (so E • co — co). Thus the vehicle 
angular momentum with respect to its mass center is 


~ c X 


j (p 4- c) dm 4- 


J p X p dm — c X 



H ~ | • co 4* Ql Zc X c 4 


j p X p dm 


(104) 


+ jp dm X (co X c) 


4- /pX(coXp) dm -1- cxjox J 

(p 4- c) dm J 

The mass-center definition requires 

(98) 

j (p 4- c) dm — 0 

(99) 

so 


J p dm ~ —cj dm — — - cQ 77 

(100) 


and 



dm 



dm == 0 


( 101 ) 


It should be noted that when c is due only to append- 
age deformations, assumed arbitrarily small, the product 
Q??c X c is of second degree and therefore negligible. For 
the vehicle illustrated in Fig. 9, the mass center of the 
appendage A 2 is nominally fixed in B,, so the large rela- 
tive motion permitted between these components does 
not contribute to c. The mass of the translating body B , 
makes the contributions of this damper body to c small 
in an engineering sense (relatively small, but not infini- 
tesimal), and for certain purposes (e.g., stability analysis 
of balanced tri-spin motion) the damper translation may 
be assumed arbitrarily small. Thus the contribution of B 2 
to QBc X c may also be negligible, For the vehicle illus- 
trated, only the rotation of A 1 on its hinge might con- 
tribute materially to c, and even this might be only an 
occasional influence. The term Q??c X c is preserved in 
this derivation for generality, but it should be recognized 
that it should be retained only in exceptional applications. 

The rotational equation (96) requires the inertial dif- 
ferentiation of the angular momentum (Eq. 104), yielding 
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(with dyadic differentiation as in Fq. 70) 




T ~ I • w 4- I • (o + )]fc X c ~i / pXp dm 


I *1 

(i + uXI'U + I • <a 4- Qli [c + 2 co Xc + oX(o)Xc) + tbXc] Xci — ^ / p X p rfm (105) 


Note that Euler’s equations emerge here for the special 
case of a rigid vehicle, since then only the first two terms 
on the right are nonzero. Another special case of interest 
is that of the nominally stationary base B 1 , for which to 
is infinitesimal and Eq. (105) linearizes to 

T = I'o 4- jp Xpdm (106) 

The challenge of Eqs. (105) and (106) is the replace- 
ment of the integrals by more convenient functions. Rela- 
tive derivatives p and p will exist only for mass elements 
that move relative to rigid body B u so the separate con- 
tributions of the rotor B ;l , the damper Bu, and the flexible 
appendages A 1 and A 2 can be calculated in turn. 

For the rotor, the relevant integral is simply 

a 1 d C 'd ' 

p X p din ” / p X (fl X p) dm — -jj (J • SI) 2 h 

(107) 

where £1 is the angular velocity of B a relative to B,, J is 
the inertia dyadic of the rotor, and h is the angular mo- 
mentum of the rotor B a relative to the base B,. 

For the transluting body Bu, the integration is trivial, 
since all parts of Bu have the same velocity p in B,. Thus 
if m is the mass of B>, and its position relative to O is 
given in the B, -.fixed vector basis b,, bu, b,„ by Bb a 4- |b , , 


with b constant and £ the single translation variable, the 
integral becomes 

>d r 0 1 d 

j pXp dm = -jj(mb£ bu) = mb£b, 4- mb£ co X bu 

(10S) 

Finally, the contributions of the appendages to this 
integral must be sought. Return for this purpose to the 
concept of the discrete parameter model of the append- 
age, and the definitions of R, r 8 , and 11 * in Eq. (36) and 
Fig. 4. In addition, let p 8 represent a generic position 
vector from the mass center P of sub-body A,, to a typi- 
cal differential mass element in A s . The appropriate inte- 
gral then becomes, for a typical appendage A (here either 
A 1 or A"), 

||pXp dm = j (R + r* + u* + p 8 ) X (r 8 4- u s 4- p 8 ) dm 

(109) 

Yet another application of the vector differentiation 
formula in Eq. (3S) is required for the interpretation of 
Eq. (109), and this can be done conveniently by writing 
the individual deformation vectors in terms of appropriate 
vector arrays (as in Eqs, 37 and 72). Recalling that is 
the angular velocity of (a) relative to {b} for a given 
appendage A, and (a} r /? 8 is the angular velocity of A* 
relative to A, one can write Eq. (109) as 


p X p dm ™ / (R + r* + u 8 + p 8 ) X [£!" X (r 8 4- u 8 ) + {a} r « 8 4- (£2“ -h {a} 7 '/? 8 ) X p s ] dm 


JA 


R 


X {a} 7 ' f u* dm + f r 8 X {a} 7 '« 8 dm 4- R X f Of 1 X [ (r 8 4 u 8 ) dm 

Ja ja L Ja 

-I- jf r 8 X [Q« X (r 8 -I- u 8 )] dm 4- t u 8 X (G« X r 8 ) dm + jf p 8 X [(&« + (a) r /? 8 ) X p 8 ] dm (110) 


since in every term the integral over A can be replaced by the sum of n integrals over the individual sub-bodies 
A t , ■ • • , A M , and 



~ 0 
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(Remember that u s and {a}' r are infinitesimal.) As for the rotor equation (107), the last term in Eq. (110) can be 
written in terms of the inertia dyadics I s of the sub-bodies. Every other integral in Eq. (110) permits factorization of 



m s 


Thus that equation becomes (with S denoting the sum over s ranging from 1 to n) 


p X p dm = R X {a } T %m s u a + Sr® X u" 

+ R X [£2 U X Sm., (r* + u®)] + tm s r s X [£2 n X (r s + u®)] + Sw»* u® X (£2° X r’) + S I* • {a} 7 ' /3 s (111) 


Equation (105) requires the inertial derivative of Eq. (Ill), which is 

f P X p dm — <o X [R X {a} T $m 8 «®] + R X [{a} r %m» i i® + (co 4- £2“) X {a} 7 ' Sm s t/®] 

+ (co + £2 a ) X Sr® X {a} 7 ’ m s u® -I- Sr® X [{a} 7 '?n s fl s + (co + £2") X {a} 7 ' m s (?] 

+ (w X R) X [£2 n X Sm s (r* + u»)] + RX [(Q“ + wX £2«) X tm s (r® + u«)] 

+ RX[Q”X Sm, (a) 7 ’ «*] + R X {£2 fl X Sm, [(co + £2«) X (r® + u®)]} 

+ Sm s [(co X £2") X r s ] X [£2 n X (r® + u*)] + Sm s r* X [(£2" -f co X £2“) X (r® -!- u®)] 

+ Sm,, r® X {£2« X [(co + £2") X (r® + u*)]} 

+ Sm,, {a} 7 ' u* X (£2" X r*) + Sm s [(co + £2") X u®] X (£2" X r®) 

+ Sm, u® X [(Q“ + co X £2") X r«] + Sm s u® X {£2« X [(co X £2") X r*]} 

+ SI® • [{a} 7 ' (S® + (co + £2") X {a} 7 ’/?®] + (to + £2°) X I® • {a} 7 ’/3® - I® X (to + £2«) • {a} 7 ’ /3® (112) 

This unwieldy expression, repeated for each appendage attached to body B u must together with Eqs. (107) and 
(108) be substituted into the rotational equation of motion (Eq, 105) in the general case of arbitrary base motion. Fur- 
thermore, the contribution of appendage deformation to I in Eq. (105) would have to be calculated. In most engineer- 
ing applications, this degree of complexity is unwarranted. In the foregoing, only appendage deformations have been 
restricted, and these have been assumed to remain infinitesimal. In fact, they are at best small in an engineering sense 
only, and deformation velocities and accelerations probably exceed, in most space vehicles, those due to prescribed 
relative motions (such as £2" and £2 a ). It is thus most reasonable in an engineering application to linearize in these 
latter variables also, although it should be understood that the resulting equations lack the rigor that can be claimed 
(however artificially) for equations such as Eq. (105), linearized in deformations only. With linearization in £2" and £2", 
Eq. (112) becomes 

J { P ^ P d™ " w X [R X {a} r Sm s «*] + BX [ {a} r ti® + co X {a } 7 'Sm,,t<®] 

+ co X (Sr® X {a } p nMi*) + Sr® [{a) 7 'm„u* + co X {a} 7, to*m®] 

+ (co X R) X (£2 n X Sm., r») + R X (£2« + co X £2«) X Sm« r® + R X [£2« X (co X r»)] 

+ Sm, r® X [(£2 n + co X £2«) X r®) + Sm„ r® X [£2« X (co X r®)] 

+ SI® • [{a} 7 ’#® + co X (a) r + co X I® • {a} 7 ’ /3* - I® X co* {a} 7 '/?® (113) 

This result, when substituted into Eq. (105), would be appropriate for the simulation of the vehicle of Fig. 9 in the 
unlikely mode of motion in which the appendages are undergoing gradual change of orientation with respect to B,, 
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while the inertial angular velocity co of B { is large. More useful forms of this result appear either when is iden- 
tically zero, which leaves 

\ si r 

pXp dm = co X [R X {a} 7 ' 2m s tV] 4- R X [{a} 7 ' 2w s u* 4- co X {a} 7 ' Sm s «*] 

4- co X (Sr s X {a } T m s u s ) 4- Sr" X [{a}' 1 ' m s u s -+ co X {a} r m g t/ s ] 

4- SI* • [{a} 7, jS s + co X {a} 7 ’ /3 s ] 4- co X I* • {a} r /?* - I* X co • {a} 7 ' /3* (114) 


or when co is also assumed to be small and included in the linearization, with the result 



(115) 

These two special cases may be of value in practical applications, since they accommodate the tri-spin vehicle with 
flexible appendages on the de-spun body, and also the dual -spin vehicle with appendages on the spinning body. Even 
further specialization may be warranted, since spinning vehicles with flexible appendages are usually simple “spinners,” 
rather than dual-spin vehicles, and a de-spun body with flexible appendages is more apt to be part of a dual-spin vehicle 
than of a tri-spin vehicle, in which case the relative motion is usually zero. 


The essential differences among the expanded forms of the rotational equation (105) are simply differences in the 
size of the equations and the amount of bookkeeping involved in writing them in matrix form. Because the purpose 
here is simply to illustrate the structure of the vehicle equations, attention is restricted henceforth to a dual-spin vehicle 
with a rigid, symmetric rotor, a nutation damper, and a single flexible appendage without rotational capability relative 
to its base (&" == 0 and C s=e E so {a} = {b}). The platform angular velocity co is unrestricted, so the equations to follow 
may easily be specialized either to the simple spinner with flexible appendages (by omitting the rotor) or to the dual- 
spin vehicle with an appendage on a de-spun body (by linearizing in co), By taking both of these steps (omitting the 
rotor and linearizing in co), one obtains equations of motion of a space vehicle with three-axis active reaction-jet control 
in an inertially stationary nominal orientation. 


By inspection of Eqs. (105), (107), (108), and (114), the required vector-dyadic equation is obtained as follows: 

T = l*co4'CoXI , co-)-l , co4-h4- mb (£ba 4~ |co X b 3 )+ u X [R X {b} 7 ’ %m s t<*] 

4 R X [{b) 7 ‘ 2m, s ff* 4- co X {b} 7 ’2m# «*] 4- co X [Sr* X {b} 7 ' «*] 4- Sr* X [{b} 7 ’m s «* 4~ co X {b} 7 'm s u*] 

4- $1* • [{b} 7 ' J'+uX (b) 7 ' jp] + co X SI* • {b} 7 ' /3 s - SI* X co - {b} 7 ' /?* (116) 


Note that the mass-center motion c has been assumed to 
be infinitesimal, in accordance with the indicated special- 
ization. 


Terms involving I in Eq, (116) are time-variable, since 
the inertia dyadic of the composite vehicle depends on 
its deformations. By definition (see Eq. 103), the vehicle 
inertia dyadic for point O is given by 

I " j (P * P E - PP) dm (117) 


Let the generic symbol p* be the value of p when the 
vehicle is undeformed, so that the inertia dyadic of the 
undeformed vehicle is 

I* = j (p* • p* E — p'" p*) dm (118) 

For all points of the vehicle except those in the flexible 
appendage and the nutation damper mass, p is identically 
p*. For points in the translating damper body Bn, 

p = p* 4- «b* = bb, 4- £b, (119) 
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(122) 


and for the mass center P s of body A s of the appendage, 

p = p* + u® = R + r s 4- u s (120) 

In addition to the contribution of the translational de- 
formation u s to the deformed vehicle inertia dyadic !, 
there is a contribution from the rotational deformation 
P® for each sub-body If the vector array {a®} consists 
of a dextral orthonormal set fixed in A s , then the inertia 
dyadic I s of A s referenced to the mass center P s may be 
written in the {a®} basis as 

I® = {a®} 7 ' 7® {a*} (121) 

If now (a®) is aligned with {b} when the appendage is 
undeformed, then the deformation {b} 7 ' /3 s yields the unit 


I = I' 6 — mb£ (b a bi 4- bib 3 ) + %m s [2 (R 4- r®) • u® E — (R -t 
Noting that the identity dyadic is 


vector relationship 

{a®} = (E - iS®) {b> 

since E — /? is the linear approximation of the direction 
cosine matrix generated by the orthogonal 1-2-3-type ro- 
tations /?,,/?-, /?»• Thus in terms of the vector basis {b}, 
the inertia dyadic for body A s becomes 

I® = {b} 7 ' (E - ft) T P (E - ft) {b} 7 ' 

= {b} 7 ’ (E + is*) 7* (E - ?*) {b) 7 ’ 

= {b} 7 ' 7* {b) -1- {b} 7 ' (ft* 7* - 7* is*) {b} (123) 

The second term is therefore the contribution of rotation 
/S® to the inertia dyadic I of the vehicle. 

From Eqs. (117) through (120) and (123), it follows that 
the vehicle inertia dyadic is 


r®) u® — u* (R 4- r*)] 4- 2 {b} 7 ’ (/?* 7® — 7® /?) (b) (124) 


E = {b} 7 ’ {b} 


(125) 


so that, for example, 


R • u* E = R T {b} • {b} 7 ‘ ti* (b} 7 ’ {b} = R 7 'ti* {b} 7 ’ {b} = {b} 7 ' R 7 'u* (b) 


one may write the inertia dyadic as 
I = I’" - {b ymb$(E l E* T 4- E*E XT ) {b} 

4- {b} 7 ' 2m* [2 (R + r*) 7 V E - (R + r») ti* 7 ' - ti* (R + r*) 7 '] {b} + {b} 7 ’ 2 (ft* 7* - 7* ft) {b} 

0 

The derivative I is readily obtained from Eq. (126) as 

f = {b} 7 ' (- mb* (E x E yr 4- E S E XT ) 4- 2m* [2 (R + r*) 7 ’ ti* E - (R + »•*) ti* 7 ’ - ti* (R + r*) 7 ’] + 2 (j§* 7* - 7* ft ) } (b) 


(126) 


(127) 


Finally, the complete and explicit vector-dyadic equation of rotational motion can be obtained by substituting Eqs. 
(126) and (127) into Eq. (116). In the process, all vectors and dyadics are written in terms of the vector array {b}, 
and the matrix representation of cross-products is employed (see Eq. 44). The result is 

[h} T T = {b} 7 ' {7*('o - mb$ (E'E* r + E*E' r ) «» 4- 2m* [2 (R 4- r*) T u* E-(fl + r>) u* - it* (R 4- r*) 7 '] & 4- 2 (ft 7* - 7* ft) <!> 
+ w7* <0 ~ 7*mb$ (E x E iT 4- E S E XT ) <o + 32m* (2 (R + r*) 7 Tt* E — (R + r s ) u sr - ti* (R + r*) 7 ’] <o 
+32 (ft 7* - I* ft*) w - vM (E'E™ + E a E 17 ') <a + 2m* [2 (R + r s ) T iV E-(R4 r») ti* 7 ' - ti* (R 4- r 8 ) 7 ’] <a 
+ S(jS 7* — 7* (8 s ) w 4- A 4- <£h 4- mb’(E“ + mb$$E 9 4- wRStn® ti* 4- R2m* it* 4' R32m* ti* + 32r* m® ti* 

+ 2? (m* it® + am* ti*) + 27* (/3* 4- 3,6*) 4- 327* /§* - 27* <sft } (128) 
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Symbols T, h, and Z* are new, being defined here by their context. The identities 


I s X co- {b} 7 ’/? 8 = I s -to X {b} T /3 s = {b} T I s Sje s 


and 


h = h + coXh = {b } T h + {b} T 7Sh 
have been employed in writing Eq. (128). 


The matrix equation isomorphic to the vector equation (128) can be written, by inspection, simply by dropping the 
vector arrays (or by pre-dot-multiplying by (b} r ). This matrix equation is more useful when written in terms of the 
total 6n by 1 deformation matrix q (defined in Eq. 83), making use of the operators % EQ and (see Eq. 85). The result 
then adopts the form 


T — P w + «Z* to + h + 3/i + 2 M (% E o R + r) T q% E0 <i» ~ * S£ 0 M (S fi0 R + r) q r Seo <I> 

- SL Mq (Sjjo R + r) r Sso « + S£ 0 (qM - Mq ) So* & - (E'E™ -f E s E vr ) + mb ($E f £ar) £ 2 

- mb ($E + |S>) (E 1 E iT + E 3 E 17 ') *> + SR Mq + R 2 7 : o Mq + RSSJ 0 Mq +ZZ r t; JMq + %lJMq 

+ Mq + SS*. Mq + S^ a rS ff «S SJ 0 Mq + Sj. 0 MSaifSS^ q - MS o*SSJ, q + qMS«s« 


— yv 

-‘OK 


s MqS ok <*> + qMSot- w — S S£ 0 qMS ok « + 2SSJ 0 M [(Seo R + r) r q] Sko « 


w M (Sko R + r) q 7 ' Sko <» — oSS 7 )„ Mq (Seo R + r) 7 ’ Seo w + 2 S£ 0 M (Seo R + r) r c'/Seo <» 

- SL [(Seo R + r) q 7 ’ - q r (Seo R + r) 7 ’j Seo «* 


(129) 


The derivation of this equation has been the principal 
objective in this section. Interpretation may be facilitated 
by restricting attention to the special case of a dual-spin 
vehicle with an appendage on a nominally de-spun plat- 
form. This case is significant from an applications view- 
point, and yet the assumption that <o is infinitesimal 
renders Eq. (129) quite interpretable. Since to is small, 
it can be replaced by the matrix 6, where 8 — [0 X 0 3 ] r . 

The result is simply the linear equation 

T = I* 9 + h — hd + mb^E' + RS£ 0 Mq 

+ 2^ 0 rMq + S 7 K Mq (130) 

Since Z * is the inertia matrix of the total vehicle as a rigid 
body, the equation T — I* 9 would suffice if the rotor were 
fixed (nonrotating) and the appendage were rigid. The 
torque applied to the body Bj in accelerating the relative 
angular rate of the rotor is — h, and the “gyroscopic stiff- 
ness” supplied by the rotor is in —hi), The term mb§E- 
reflects the “inertial torque” due to damper mass acceler- 
ation, and the last three terms reflect the effects of ap- 
pendage deformation, 


The more general Eq. (129) has the major added term 
«/* ( 0 , which introduces the “Euler coupling” or “gyro- 
scopic coupling” of the vehicle as a spinning rigid body. 
In addition, Eq. (129) includes a multitude of nonlinear 
terms that reflect the influence of the vehicle deform- 
ability. They can be categorized as various kinds of accel- 
eration terms (centripetal, coriolis, etc.), but perhaps they 
are best understood by reviewing their origins in the 
vector-dyadic equation (116). 

The matrix equation of general rotation of the vehicle, 
Eq. (129), is compatible in its assumptions with Eq. (95), 
which defines the appendage response to a given arbitrary 
base motion. These equations do not yet constitute a com- 
plete set, however, since they include (in addition to 
unspecified external forces and torques) the unknown 
motion variables X , $, and h. 

_ * * 

The 3 by 1 matrix X is the matrix in an inertial vector 
basis of the inertial acceleration vector of the vehicle 
mass center. This vector is easily related to the resultant 
applied force F by applying Newtons second law, i.e,, 

F = QftX = QR {i} r X 
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If F is written in terms of the body-fixed vector basis {b}, 
related to {i} as in Eq. (53), this equation becomes 

{b} r F = 111 {i} T X = m {b} T ex 

or in matrix terms 

F = (131) 

Newtons second law also provides the equation of 
motion of the damped linear oscillation representing the 


nutation damper, except that only the b r component of 
the vector equation of translation is required. The vector 
result is 

M (X + c + bi> 3 + hi + fbi) • b ! = F Si • b x = -fcf - 4 

(132) 

where k and d are scalar constants of the spring and dash- 
pot, respectively. The mass-center displacement vector c 
is available from Eq. (60) as 


c -{b } T ZnMq - (m4/Ml)b x = -{b) r [SsaMq + m$E'/Ml] 

so that 


•{b} r [Sbo p.(j + m£E\ Ml 4- 25" (2i,- 0 4 4- m%E 1 QH) 4- w (2g„ nq 4 m$E 1 /Ml) 4- <35? (Sbo pq + m$E l /Ml)] 


(133) 


Similar operations on b„ — {b} 7 ' E :i and b t = (b} r £ l permit the expansion of the term in parentheses in Eq. (132), and 
the b v dot-product then provides the scalar equation 

in (1 — rn Ml) £ 4 (d- — 2mE lT wE 1 m Ml) £ 4 [k — mE lT (Z 4- <5S) (1 — rn/Ml) E l ] £ = 

mE vr [ — C-)X — b (<L 4- 5£>) E s 4- p.q 4- 2 S2 £0 ^q 4- (i 4- SS) Ssu jito] (134) 


It should be noted that the ratio in Ml is too small to 
’’.'arrant retention when B t is a nutation damper, although 
for other applications of translating second bodies this 
may not be the case. 

Finally, an equation of motion must be provided to 
determine the remaining variable 7i, the angular momen- 
tum matrix of the rotor B ;t relative to body B u expressed 
in the vector basis of B t . Because the rotor spins about 
its axis of symmetry, its relative angular momentum can 
always be expressed as the product of a scalar 7<£ and 
a unit vector fixed in B u where *7 is the moment of 
inertia of the rotor about its symmetry axis and <+> is the 
spin rate of the rotor B 3 relative to the base B x . The 
orientation of the rotor in B v is here selected so that the 
rotor axis parallels b.j, merely for analytical convenience. 
Thus the relative angular momentum h may be expressed 
as 

h = 7<£E 3 (135) 

Thus a single scalar equation must be derived for the 
determination of the scalar unknown <£. This is readily 
obtained as the scalar Euler equation corresponding to 
the axis of symmetry, namely, 

r = 7 (<j> 4* « 3 ) = .7 (<£ 4- E 3T w) (136) 


where r is the magnitude of the torque applied to the 
rotor along its axis of symmetry. In general r is some 
combination of bearing friction and the motor torque 
established by a control law. Complete simulation would 
of course require an auxiliary equation for that control 
law. and except for very preliminary calculations the 
control law would be nonlinear. 

In summary, Eqs. (136), (134), (131), (129), and, (95) 
constitute a complete set of dynamic equations for the 
system. The unknowns are included in </>, £, X, (-), o>, and q. 
(Equation 135 must be substituted into Eq. 129 before 
collecting these equations, in order to eliminate h.) To 
complete the simulation, a set of kinematic equations 
relating (-> and a must be included. The particular equa- 
tions used will depend on the attitude parameters chosen 
to define the orientation of B x in inertial space (e.g., direc- 
tion cosines, Euler parameters, attitude angles, etc.). 
These equations are available from the fundamental kine- 
matic relationship (Ref. 16) 

~ = 0 0 r = _ 0 (137) 

but more suitable forms of the kinematic equations can 
be obtained for specific parameter choices (e.g., see 
Ref. 23). 
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It is then conceptually possible to collect the five dy- 
namic equations (136), (134), (131), (129), and (95) to- 
gether with kinematic equations from Eq. (137) and 
control laws for the rotor torque r and the vehicle torque 
T, and with this set, to accomplish the simulation of a 
dual-spin vehicle with a flexible appendage and a nuta- 
tion damper attached to a base undergoing arbitrary 
motion. By specializing these equations, one could simu- 
late a simple spin-stabilized vehicle, or a vehicle with 
three-axis reaction-jet control. By generalizing these equa- 
tions, proceeding from vector-dyadic equations developed 
here in detail, one could with straightforward labor simi- 
larly simulate a vehicle controlled by several rotors 
(“momentum wheels”) and having articulated elastic ap- 
pendages capable of substantial relative motion. All of 
this is, however, only conceptually possible. The complete 
set of equations is available, but if they are to be applied 
to a real space vehicle represented by a reasonably accu- 
rate mathematical model, the dimension of the resulting 
equations is prohibitively high, even for machine com- 
putation. Remember that Eq. (bo) alone consists, for 
practical problems, of hundreds of second-order scalar 
differential equations, and that the entire system of equa- 
tions is nonlinear and highly coupled. If these equations 
are to have practical value, they must be subjected to 
coordinate transformations that accomplish substantial 
uncoupling of the equations and permit the analyst to 
exercise judgment in restricting his attention to a subset 
of the present set of perhaps hundreds of coordinates, in 
order to justify working with a more reasonable number 
of differential equations in the simulation of the vehicle. 

D. Coordinate Transformations 

1. Transformation rationale. The objective of this sec- 
tion is to investigate the possibility of simplifying the 
equations of motion by introducing linear transformations 
for some of the variables. The transformations are applied 
to the appendage deformation matrix q only, leaving the 
discrete coordinates of the base B u the damper B 2 , and 
the rotor B a unchanged, The more remote possibility of 
finding a transformation that may fruitfully be applied to 
the entire system of coordinates is treated in Section IV. 

Of course, an infinite variety of coordinate transforma- 
tions may be devised, even within the linear class (to 
which attention is presently restricted). Any transforma- 
tion that provides a one-to-one mapping from one coordi- 
nate system to the next yields equations that could 
theoretically provide a correct simulation, but obviously 
not all such transformations are advantageous. It may be, 
however, that more than one useful transformation will 
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be available, and then the appropriate choice may depend 
on individual standards of utility. 

The basic difficulty in using the equations of motion 
derived earlier (Eqs. 95, 129, 131, 134, and 136 or their 
equivalent) is simply their dimension. A certain amount 
of nonlinearity and coupling of equations may be un- 
avoidable. Therefore the primary test of the utility of a 
given coordinate transformation is the degree to which 
it permits the truncation of the coordinate matrix, and 
the consequent reduction in the dimension of the system 
of equations. 

Coordinate truncation is never a completely rigorous 
process, since it results in an incomplete and imprecise 
indication of the response of the mathematical model to 
its dynamic environment. As a practical matter, however, 
it must be recognized that the mathematical model is but 
an imperfect approximation of the vehicle being studied, 
so it is sophistry to argue that coordinate truncation nec- 
essarily degrades the simulation of a real vehicle. Some 
degree of truncation may be appropriate even when hie 
most realistic efficient simulation is sought. 

There are no entirely satisfactory analytical procedures 
for determining the degree to which a given matrix, of 
coordinates can safely be truncated, nor can the relative 
acceptability of truncation of alternative coordinate sys- 
tems be unequivocally established. It should be noted, 
however, that when a coordinate transformation pro- 
vides a completely uncoupled system of scalar differential 
equations, one can obtain the solution of these eq' ations 
precisely for each scalar coordinate in turn, without con- 
sidering the influence of other coordinates. Therefore one 
can truncate the coordinate matrix without sacrificing 
the rigorous validity of the solutions for the coordinates 
retained. The determination of the number of such coor- 
dinates to retain for acceptable simulation is still a matter 
of engineering judgment, but the validity of the solution 
for those retained is related directly to the degree to 
which the transformed equations are uncoupled, Accord- 
ingly, the desirability of a given coordinate transforma- 
tion is measured here in terms of the degree to which it 
uncouples the system of differential equations employed 
in the vehicle simulation. 

The total system of equations required for the vehicle 
simulation includes some equations in which the append- 
age ch'foTmatlon terms are of central importance (e.g,, the 
appendage equation 95), some equations in which terms 
in q are generally of secondary importance (e.g., the nuta- 
tion damper equation 134), in which q is always multi- 
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plied by the small constant matrix fi), and some equations 
in which q does not appear at all (e.g., the rotor equa- 
tion 136), In the quest for a transformation that to some 
degree uncouples the equations, it is therefore reasonable 
to concentrate on a subset of the system of differential 
equations, excluding from consideration the equations of 
motion of rotor and damper, the kinematic equations 
relating o> and 0, and the control system equations. (Note 
that the last of these will include q when sensors or 
thrustors are mounted on flexible appendages.) This deci- 
sion to concentrate on certain equations clearly marks a 
compromise with the objective of securing uncoupled 
equations, and, in fact, this objective became unrealizable 
with the announced restriction in this section to transfor- 
mations of the appendage deformation coordinates q only. 
The Anal equations must certainly remain coupled by the 


discrete coordinates of B 1} B>, and B ; „ and by ignoring 
certain equations while transformations are sought, one 
accepts as well the likelihood that these equations will, 
after transformation, remain coupled in the new deforma- 
tion variables. 

Attention is now restricted to Eqs. (95), (129), and (131), 
which are respectively the appendage equation, the ve- 
hicle rotation equation, and the vehicle translation equa- 
tion, The last of these is trivial, permitting the solution 

eX = F/Qn (138) 

This substitution into Eq. (95) with the substitution 
e = —m£E 1 /(y)'i to accommodate the damper, provides 
the revised appendage equation 


M (E - v, ;n ll n M/m ) q + 2 (M [(S*o «>)"- S*, & M/m ] + M (S„« u>)~ + (So* ^M- u>)~+ D') q 

~ IM [(2/jn w)~— <0 Z r m M/m) + M [(2/ JO «‘))~(Sho«>)~ — Smoww M/m] + K) q ~ 


■M (X 0 /.’ ~ Xi:n R ~~ M [Sen omR — (S« n«») (Xnn '<>) >’] — * (S<//j «) M (Sn/j ">) 

~ MS /jo F/m + MS/jo m‘£E x /m + X 


(139) 


This equation consists of 6n second-order scalar equations and can be written as a matrix equation of the structure of 
Eq. (84), or somewhat more explicitly as 


M'q + D'q + G' q + K' q + A' q — V - -M(S os 


where N° represents the nonlinear terms in «> due to cen- 
trifugal forces on the appendage, Here G' depends lin- 
early on < 0 , A' depends linearly on d>, and K' depends non- 
linearly on «> (involving only second-degree terms). As 
noted in the discussion of Eq, (84), M', D', and K! are 
symmetric and G' and A' are skew-symmetric, 

The remaining equation of primary interest, Eq. (129), 
consists of only three scalar equations, although many of 
its terms involve matrix operators that accomplish sum- 
mations over On elements, The appearance of Urn variable 
q and its derivatives in Eq, (129) is again linear, and again 
ilt appears only in the first power and «i in the second, so 
it would be possible to write this equation in the form 
of Eq, (84). This would be misleading, however, since the 
coefficient matrices would be rectangular (3 by On), and 
the role of Eq, (129) as the vehicle rotational equation 
would be obscured. To emphasize the primary physical 
significance of this equation, it is rewritten in symbolic 
form as 


S/joH - rS 0 #)A + N' ~ MXm F/m + M% m m$E'Wl + X 

(WO) 


Id) + dllli) + (/ -|~ G — h) < I > — 

-(Sj # + R Sj # + ZJJ) Mq - h - rnb$E" + T 


(141) 


This form of the vehicle rotational equation can be con- 
structed from Eq, (129) quite mechanically by defining 1 
as the coefficient matrix of d>, observing that this same 
matrix permits whu to accommodate all second-degree 
terms in «>, and then defining G to include ajj terms in the 
coefficient matrix of w not absorbed by I or h, Much more 
insight into the significance of 1 and G may be gained by 
inspecting Eq, (116), which is the vector-dyadic antece- 
dent of Eq, (129), With this examination comes the real- 
ization that J is the time-variable inertia matrix of the 
total vehicle in basis {b}, referred to the vehicle mass 
center, The term Gm accommodates the "inertial torques” 
applied to the vehicle due to the coriolis accelerations of 
appendage and damper, so it has the same physical ori- 
gins as the matrix G 1 of Eq, (140), It should be noted that 
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the coefficient of q in Eq. (141) is the transpose of the 
coefficient of <» in Eq, (140), i,e„ 


[(*?;, i i ^„r) m r - a 

■ M (v nw - v wn R - 7v W(| ) 

(142) 


since a symmetric matrix is its own transpose and a skew- 
symmetric matrix is the negative of its transpose. 


Equations (.140) and (141) are coupled with each other, 
and Eq, (141) is also coupled with the equations of motion 
of rotor and damper, Nonetheless, physically based argu- 
ments can be introduced that tend to justify the use of 
a transformation on q, which uncouples the new deforma- 
tion variables in Eq, (140) only, leaving linear combina- 
tions of these variables in each of the scalar equations of 
Eq, (141), Truncation of the new deformation coordinate 
matrix would then permit substantial reduction of the 
number of scalar equations accepted from Eq, (140), and 
this would provide exact solutions for the coordinates 
retained only insofar as m (as obtained from Eq, 141) is 
uninfluenced by the excluded coordinates, This seems 
quite an acceptable approximation, since in general the 
appendage vibrations act only as perturbations on the 
motion of the controlled vehicle, which is dominated by 
the control torque T in Eq, (141). 


Finally, then, attention focuses on the quest for a coor- 
dinate transformation that uncouples the deformation 
coordinates in the On equations of Eq, (140), Only the 
homogeneous equation 


No general procedures are available for the transfor- 
mation of Eq. (143) into a system of uncoupled equations, 
This can be accomplished in general terms only when the 
coefficient matrices are constant, which requires that the 
nominal value's) be constant, This is surely the most com- 
mon situation of interest in spacecraft control, but it 
excludes certain appendage deployment and transient 
operations problems, 


2. State-equation modal analysis, Consider then the 
special case of Eq, (143) for which A/', D\ G', and K! are 
constant, implying AheO, This second-order matrix equa- 
tion can always be represented by a first-order equation 
(the state equation), This is accomplished by defining the 
12n by 1 matrix 


and writing 
where 



Q - BQ 


(144) 


(145) 


11 vU 


0 


L -( A /')- 1 K ' 


E 

(a ry (g~~d')~ 


(146) 


Here B is a 12n by 12/j matrix, and each of the partitions 
shown in Eq, (146) is of dimension On by On. Equa- 
tion (145) has exactly the same content as Eq, (143), since 
it contains this equation in addition to the identity q " q. 


Linear, constant-coefficient differential equations such 
as Eq, (145) always have solutions of the form 

Q^<|>m e W (147) 


M'q i I)' q I (V q 4 K' q f /V q = 0 (143) 


is relevant to this search, 

Although matrices A/' and D' are constant, matrices 
O', A', and K' may vary with time, since they depend on 
w and m, Of course <» and <’,• are strictly unknown, but it is 
quite reasonable in the present context of space vehicle 
attitude-control-system simulation to assume that m varies 
only slightly from the nominal value (say ,To(t)) which is 
the objective of the control system, Under this assump- 
tion, formal linearization removes the unknown <» 
from the homogeneous Eq, (143), since this small quantity 
appears in that equation only as a factor of q or q, Thus 
one may argue quite generally that Eq, (143) is a system 
of linear second-order equations with coefficients depend- 
ing explicitly on time, 


where \ m is a scalar and 4> HI is a 12a by 1 matrix, which 
from Eq, (144) must partition into a On by 1 matrix <j> m 
above a On by 1 matrix <j> m A m , i.e,, 


#* = [*? (tjtf v) ■ ■ ■ m.uv 


~ </,"• - 

(148) 


The validity of the solution indicated in Eq, (147) can 
be confirmed by substitution into Eq. (145), yielding 


a„, <i> m - m m 

since 0, This result may be written in the conven- 
tional form of an eigenvalue problem, 

(B - km E) <F" « 0 (140) 
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so that Am is an eigenvalue and <I»' H an eigenvector of B. 
Tlie existence of a solution as in Eq. (147) follows from 
the existence of at least one eigenvalue of B and its cor- 
responding eigenvector. 

In scalar terms, Eq. (149) represents a set of 12n homo- 
geneous algebraic equations in the 12n + 1 unknowns 
‘I - * • , «!>!»„ and A,„. Nontrivial (nonzero) solutions for 
<I> m exist only when A m takes on certain characteristic val- 
ues or eigenvalues such that the determinant of coeffi- 
cients is zero, by Cramer's rule. In other words, values 
of A m must be selected that satisfy 

| B - Am E | r 0 (150) 

There arc 1.2 /j solutions A t , • • • , Aia« for Eq, (150), as 
may be confirmed by expanding the determinant into the 
polynomial 

b VJn A 1 "' 1 H ft,..,, , A ia ” 1 -I T b. A" + ft, A T b, ■■ 

(A — A,) (A - Aa) • ’ • (A — Amu) r 0 

(151) 

The matrix B and the related cofficients b lh • • • , l) v , n are 
real numbers for equations that stem from the dynamic 
equations of Eq. (143), Therefore the complex eigenvalues 
among A,, 1 ■ • ,A, B « occur in complex conjugate pairs, 
Subscripts may be assigned so that Am,.,,,, ~ A* (asterisk 
here denotes complex conjugate), 

Corresponding to each distinct eigenvalue there 
exists an eigenvector <l* f " that can be determined .o within 
a multiplicative constant by solving Eq. (140 A com- 
pletely unique solution is not possible, sine the sys- 
tem of equations in Eq. (149) includes on , 12«. — 1 
independent algebraic equations in the 12n . nknowns 
* * • *<!*?{,«* with Am specified by Eq, (150). \v\on the 
eigenvalues are distinct, the eigenvectors arc independent 
(see Ref, 33, pp, 184 IT,), but when there are repeated 
eigenvalues, the corresponding eigenvectors may not be 
independent. 

The validity of Eq. (149) for a particular A», and <l> m 
guarantees as well the validity of 

(Fa;£)-I'^-Q (152) 

since the complex conjugate of zero Is zero, and the con- 
jugate of a product Is the product of conjugates. Because 
K - Afiii+wi the eigenvectors corresponding to complex 
pairs of eigenvalues are, from Eq. (.152), also complex 
conjugate pairs, 


Although the solution to Eq. (145) has readily been 
found, and certain of the properties of the solution have 
been noted, these steps have not yet accomplished the 
main objective, which is the discovery of a transformation 
matrix that uncouples the coordinates in Eq. (14.3), It is 
in most cases not essential that the result he a second- 
order differential equation; it may suffice to obtain a first- 
order equation such as Eq. (145), but with a diagonal 
coefficient matrix, Although it is possible to generate a 
transformation from Eq. (143) to an uncoupled system of 
second-order equations when A/' and K' are symmetric 
and O' is skew-symmetric, it is possible to diagonalize the 
state equation (145) even more generally, The latter trans- 
formation is considered first, 


If <l> is defined as the 12» by 12n matrix whose columns 
arc the eigenvectors <!>', • * * , of the matrix B, then 
by inspection of Eq, (148) one may write 


B<|>; B ■ • 1 'I*' 8 "] 
[Ai ‘I» l ’ ' * Aisn'k 18 "] 


[«}»> • ■ • ({>'««] 


A, 


0 


A a 


o - 

Aian, 



0 " 
Aiuii _ 



It follows that, if 4» 1 exists, promultiplication provides 
the diagonal matrix 


<I> ' /M> 


"A i 

0 


0 

Aia ti 


(154) 


The existence of <l» 1 is assured if the determinant of <!' 
is nonzero, i.e,, if the eigenvectors 4' ', 1 ■ • ,«!» ,8N are inde- 
pendent, It lias been noted that this independence is 
guaranteed when Ai, * • • , Am» arc distinct, For the dy- 
namic system of Eq, (143), the independence of the eigen- 
vectors can be assured even without this restriction when 
A/\ G\ and K' are constant and A ' and D' are zero, The 
necessary supporting argument is briefly summarized in 
the following paragraphs. 
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When the eigenvectors are independent, the general 
solution to Eq, (1*15) is the linear combination of solutions 
of the structure of Eq, (.1*17), i.e., 

() V /\„.<|.»» <>•'»' (155) 

id I 

The scalar constants A m are established by initial condi- 
tions, It follows from Eq, (155) that if the real parts of 
A,, • ■ * , A.,-,, are nonpositive, the response to initial con- 
ditions providing sufficiently small values of A,, • • • , A V j„ 
is a solution Q with an arbitrarily small norm, In other 
words, when the real parts of the eigenvalues of B are 
nonpositive, the null solution of Eq, (1*15) is Liapunov 
stable, This is true even if there are repeated eigenvalues, 
as long as the eigenvectors are independent, 

If on the other hand there were repeated eigenvalues 
(say, A; A,,)) with dependent eigenvectors, the general 
solution would become 

() = A i 'I* 1 c A| ' 4 ■ • ■ I Aj 

t Am «l» ; tc*' 1 '! • • * 1 A,a„ 'I* 1 '" c w (156) 


Because the matrix C ' in Eq, (143) is skew-symmetric, 
and the symmetric matrix D ' has been deleted, there 
can exist no eigenvalues with negative parts, This is 
a consequence mathematically of the absence, of odd- 
powered terms in the characteristic equation (151), In 
physical terms, this is a reflection of the fact that the 
appendage as presently idealized has no damping, and 
consequently no attenuating solutions are admissible, 
Thus Liapunov stability of the null solution of Eq, (143) 
must mean that all eigenvalues A,, • • • , An..,, have zero 
real parts, under present restrictive assumptions, 

Consequently it can be concluded that, for this class 
of system, the eigenvectors <l>', • • • , <T»' 8 " are indepen- 
dent, and the matrix 4' 1 always exists, This means 
that the operation of Eq, (154) is always possible, and 
matrix B can always be diagonalized for the problem 
under consideration, 

In application to the differential equations of the space 
vehicle, the appropriate transformation procedure is 
straightforward, Each of the equations involving q (Eqs. 
140, 141, and 134, for example) must be rewritten in terms 
of Q as defined in Eq, (144), Thus Eq, (140) becomes 


If the real part of A, is zero, the function /r v is un- 
bounded, and the null solution of Eq, (145) is unstable, 
Liapunov stability is possible in the presence of depend- 
ent eigenvectors only if the real part of the corresponding 
eigenvalues is negative, 

Now it will be argued that the null solution of Eq, (143), 
restricted to constant A/', G\ K\ and to zero values of 
A' and /)' (and hence of the corresponding Eq, 145) is 
Liapunov stable, and that the eigenvalues Ai, • - 1 , A ly „ 
have zero real parts. These conditions have been shown 
to be incompatible with the presence of dependent 
eigenvectors, 

The function q f M'q i q 1 K'q is a Liapunov function 
(see Ref, 34, Oh, l) when A/' and K' are positive definite, 
so under these conditions the null solution of Eq. (143) 
(and hence of Eq, 145) is Liapunov stable, By definition 
(Eq, 37), A/' is positive definite. The* matrix K\ defined 
by Eq, (91), consists, wfuv w is constant, of the positive- 
definite stiffness matrix a plus a matrix provided by the 
centrifugal "force" field induced by vehicle rotation, Un- 
less centrifugal forces are so great as to exceed the elastic 
forces of the structure (in which case the appendage 
would pull apart), the matrix K' is positive definite, Thus 
under all reasonable circumstances the null solutions of 
Eqs, (1.43) and (145) are Liapunov stable, for the noted 
restrictions, 


Q BQ H L (157) 


where, in terms of 6n by 1 submatrices 0 and V (see 
Eq, 140), the 12n by 1 matrix h is [0 V] r , Equation (141) 
may be written as 


1 m 4 ~ mJm + (I + G —• /i) (i) ™ 

{0|-(sj, + 8ss. + SE. 7) M] Q - ■ ft 


EM T 
(158) 


« 

where the 3 by 12m coefficient matrix of Q is shown in 
two 3 by 8m partitions, Recall that in writing Eq, (141) 
from Eq, (129), the matrices / and G were noted to de- 
pend on q and q, so these matrices in Eq, (.158) require 
reformulation in terms of (), Similarly, the damper equa- 
tion (134) and any others involving q arc easily rewritten 
in terms of (), In practice, one would normally carry the 
transformation into first-order equations all the way, re- 
placing, for example, the second-order scalar damper 
equation by two first-order equations by introducing a 
6 by 1 matrix for U i]' r ns was done in Eq, (144), and 
similarly replacing | in Eq, (158), 


The transformation 

Q '-■■■ «l»Y (159) 


Is next Introduced wherever Q appears in any of the 
equations of the vehicle or control system, and Eq, (157) 
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is then premultiplied by <T r5 , The transformed appendage 
equation then becomes (by virtue of Eq. 154) 


y 


y + ‘i’"’ l 


A, 0 

, 0 Acjii J 

and the vehicle rotational equation (15S) becomes 


(160) 


/<!) "1 mini A (I A G •— /)) <i> “ 

[o| - (s;„ + ns;,, + sj,?) m] *r - A - mbH& + r 

(161) 


The scalar equations in Eq. (160) are not directly cou- 
pled by the deformation coordinates in Y. It should be 
noted, however, that the matrix L in Eq, (160) is deter- 
mined by the base rotational motion obtained from 
Eq, (.161), and in general, all of the scalar variables in Y 
participate in Eq. (161), Thus when the total system of 
differential equations is considered, there is still coupling 
among the deformation coordinates in Y, Truncation of 
the 12n by 1 matrix Y and the corresponding 12m by 12m 
transformation matrix <l> must therefore be based on the 
considered judgment that the coordinate coupling is weak 
and physically indirect, Such arguments support the prac- 
tice of determining coordinate truncation principally from 
the uncoupled appendage deformation equations in 
Eq, (160). It may be necessary to consider also certain 
characteristics of the control system, such as the range of 
frequencies of response, sensor data frequency filtering, 
and the location of sensors and thrusters, Discussion of 
these questions of control system interaction is postponed 
until Section V, 


The homogeneous solution for Q in Eq, (155) is, of 
course, the same as that obtained by solving Eq, (160) for 
Y when L is zero and substituting into the transformation 
of Eq, (159), Truncation of Y is normally accomplished 
by inspection of this homogeneous solution, although it 
may also be necessary to include any coordinates in Y 
that may introduce resonances due to coincidence of the 
imaginary part of the corresponding A with a driving fre- 
quency in the forcing term h in Eq, (160), Since trunca- 
tion is determined from physical considerations, it is 
important that the physical signfieanee of the variables 
Yi, > • * ,Y m be understood, This may seem difficult, 
since complex eigenvectors in A make the coordinates 
Y„ * * ■ , Tib,, generally complex, even when Q is by defi- 
nition real (see Eq, 159), The indicated homogeneous 
solution of Eq, (160) has the elements 

Y m — A m e ,w (162) 


Thus the solution shown in Eq, (155) may be written as 

Q 2 Y,„ <l» m (163) 

m i 

in conformance with Eq. (159), Although individual sca- 
lars Y,„ and column matrices A" 1 are complex, the reality 
of Q is assured by the appearance of eigenvalues and 
eigenvectors in complex pairs, Consider, for example, a 
solution for Y with initial conditions selected so that all 
values of A,„ in Eq, (162) are zero except A h and A,,,,.*, 
Then Q is 

Q” y K <[,*• + Y (I(I+a - A"" +a ' 

* ■ Ax A a T- An,,,/,- A" ,ltA ' e x "** ht 
~ Ah «l > K cA* 1 -I A,,,, c Xl '' t 

Let A + and AY 'I'' 1 I iV K and substitute 
r u ' - c 011 ' 1 (cos a/, t I /sin cr A /), etc, Then Q becomes 

- c n «> {(Aa- I A,i„+a ) (A* cos <y a I - - V h sin <r A I) 

A i (Aa ~ A,|„* a ) (r A cos ir K t A A A sin cr A /)} 

Since Q is real, A ( „ i+ a- is A* , Let Ak ss (Cl* — /£«■)/ 2, so 
that Q may bo rewritten in wholly real terms as 


0 = c nKt {C a- (A* cos <r K t — V K sin a*/) 

A D k (a* sin a h t A V K cos <)■*/) } (164) 

Thus the general homogeneous solution for Q appears in 
real terms as 


n,i 

Q - 2 {Cm (i'" 1 COS <r,„t — r m .sin <r m t) 
m-T 

A D„ i (A m sin (Tint A r« cos «t„,/)) (165) 

The 12 m scalars Ci, * • * , C,,,, and D t) • • ■ , D„„ are estab- 
lished by the initial conditions on {A For an undamped 
system, as in Eq. (143) when DY .0, the scalars m, • * • , (v,m 
are all zero, 


Now coordinate truncation can be based on the same 
rationale traditionally employed by structural dynnmicists 
in application to second-order equations of motion, Each 
conjugate pair of coordinates (say, Y m and Y mm ~ Y*,) 
corresponds to appendage vibrations of the entire struc- 
ture at a given frequency (here «,,,). In general, these 
vibrations are not in phase throughout the structure, as 
would he the case for classic normal modes of vibration 
about a rest configuration, The classic interpretation of 
the eigenvector A'" of Eq, (148) as the mode shape of the 
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vibration at frequency v m is modified by the complex 
nature of <l» m — •! iv m . As illustrated in real terms in 

Eq, (164), this is reflected as two different mode shapes, 
corresponding to two vibrations at frequency a m> with a 
90-deg phase lag, Furthermore, these mode shapes de- 
pend on the initial conditions. The coordinates in V are 
still called modal coordinates in this report, despite these 
deviations from the classic interpretation, 

The usual decision of the structural dynamicist is to 
retain in truncation the responses at the lowest natural 
frequencies, and also any others at which 'expected driv- 
ing frequencies are close to natural frequencies (with 
consequent danger of resonance), The justification for 
concentrating on lower frequency responses is often 
based on the fact that storing a given amount of energy 
as potential energy or strain energy in tho deformation 
of an elastic structure into a low-frequency mode shape 
generally results in larger displacements and base reac- 
tions than are present when the same energy is stored in 
a high-frequency mode shape deformation. To this it may 
be added that the present clastic model is unrealistic in 
its neglect of damping, and that, in fact, the higher fre- 
quency vibrations will have more damping than those at 
low frequency, This again supports tho practice of ignor- 
ing the high-frequency response (since it rapidly attenu- 
ates in the transient solution anyway), 

The appendage deformation matrix Y in Eqs. (.160) and 
(.161) is therefore replaced by the 2 N by 1 truncated 
matrix 

V: ..[Y t • ■ • YvY* ■ • ’ Y* v j T (166) 

where N is the number of modes to be preserved in the 
simulation, The transformation matrix «l> is accordingly 
truncated to the 12 n by 2 N matrix <T>, where 

fejy • - • <|>v (167) 

The equation of motion of the appendage (Eq, 160) now 
becomes 


A, 


0 




A? 


Y -I- L (168) 


0 


Am 


This matrix equation consists of 2N scalar equations in 
the complex variables Y 1t * * ' , Y. v , Y* , • * * *Y$, In com- 


putation, only the first N equations need be considered, 
since the equations are uncoupled and the second set of N 
variables is available as the complex conjugates of the 
first set, 

The rotational equation of the vehicle (Eq. 161) then 
becomes 


Jm 'f John + (/ + G — /() tit — 

[0! - (SJ K *b + 3j) (1 7) M) $7 



The matrix <7» is rectangular, and the symbol <1» ’ in 
Eq, (168) represents a left inverse of <1*. This is available 
(Ref, 35) as 


<r> = (<T> T <T>b l (!>'«’ (170) 

as may be confirmed by the multiplication 

[(‘T* 7 ’ <7>) 1 <T» r ] <T> (7> r T*) 1 (<T‘ r <T>) — e 

The inversion actually performed in Eq, (170) to obtain 
<S> 1 is of a matrix of dimension 2 N, Computations aro 
thus not as imposing a task as may have been implied by 
the presence of the 12n by I2u matrix inverse <l> 1 in 
Eq, (160), In practice, the numerical value of «(> is never 
computed, since only <T is required, 


Incorporation of energy dissipation into the appendage 
mathematical model is usually deferred, pending trans- 
formation to modal coordinates, The scalar (v A in Eq, (164) 
would be zero for the undamped system of Eq, (.143) with 
P' ~ 0, since G' is skew-symmetric, Conventional practice 
in structural dynamics is to assign values to <x u * 1 1 ,«, v 
(where JV modes are retained in the modal-coordinate 
truncation), with preliminary numerical values based on 
experience with similar structures and final values based 
on prototype tests, This amounts to rather arbitrarily 
adding a viscous damping term P' q to Eq, (143) after 
completing the modal analysis, and assuming that the 
structure of the matrix P' is such that the eigenvectors 
</)', • ■ ’ are undisturbed by this addition, Numeri- 
cal values of the elements in D' are never considered, 
since they are determined implicitly by the choice of 
m. ' ' ' »fw< (Actually, the conventional pattern is the 
assignment of a value to *= <r m rather than to <v m 
Itself, Here is called the "percentage of critical damp- 
ing," For typical spacecraft appendage structures, is 
well below 5% for the lower frequency responses,) 
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(171) 


The mechanical process of incorporating clamping in the 
appendage equations is straightforward, Equation (168) 
includes a diagonal matrix of eigenvalues A,, ■ * • , A.v, 
A*. • ’ • .A*. As calculated from the matrix B of Kq, (148) 
with jy set equal to zero, these eigenvalues will he purely 
imaginary, e,g„ Am r fo m . To incorporate structural damp- 
ing, one instead substitutes A,„ £ m < h„ 4 into Eq, (168), 

where Cm is chosen as noted in the preceding paragraph, 


It should be noted that the procedure just described 
for transformation to modal coordinates and subsequent 
truncation did not require any special properties of the 
constant matrix B in Eq. (145), except in proving the 
existence of <!> 1 in Eq, (.154). The nonsingularity of <l> 
was established by relying upon the skew-symmetry of 
matrix O' and the absence of the viscous damping ma- 
trix D' in Eq. (143), This property was used to prove that 
even in the special case of repeated eigenvalues, the 
independence of tho eigenvectors is assured, so the solu- 
tion to Eq, (145) has the form of Eq, (155) rather than that 
of Eq, (156), The proof would proceed equally well if G ' 
were zero, but when damping is present the formal argu- 
ment underlying the previous discussion breaks down, 
Although it seems reasonable that continuity would 
require the preservation of the nonsingularity of <l> in the 
presence of small damping forces (with the correspond- 
ing V <) included in Eq, 143, with O' symmetric), proof of 
this condition seems lacking, 


It can be asserted, then, that only when «(> 1 exists is 
tho procedure applied here to transform the state equa- 
tion (145) into the diagonal form applicable, It is per- 
haps worth repeating that the existence of ’ is assured 
when the eigenvalues of B are distinct, or when the 
matrix £>' is zero, This method has thus been shown to be 
applicable to any linear dynamic system with distinct 
natural frequencies, and to any undamped linear dynamic 
system, Applicability to damped linear dynamic systems 
with repeated natural frequencies seems probable, but has 
not been formally established, 


it/" Q + K" Q ~ L" 

where 


A/" 


D‘ [ A/' 

„ ° , 


K" E~ 





77 




Foss demonstrates that, for symmetric D\ the eigenvec- 
tors and *|A corresponding to different eigenvalues 
(A w A/,) have the orthogonality property 

<|,'" T A/" <|A 0 ) 

> (173) 

<!•"< ' K" <!>'• • 0 1 


Thus for suitable eigenvector normalization, diagonaliza- 
tion of Eq. (17.1) with the transformation (see Eq, .159) 



Q - <1*V 


is computationally easier than the diagonalizatien y)f 
Eq, (145), Whereas the latter procedure requires tho cal- 
culation of <l> 3 * 1 * (or <T> 1 after truncation), the former re- 
quires only the construction of <t* r (or <T> r ) to obtain the 
same result. This transformation of Eq. (171), followed 
by premultiplication by <l» 7 ', provides the same results as 
Eq, (168), namely, 


¥' M" <t>Y l <|> 7 ' K" 4>V " <l» 7 ' L" 


3, Modal analysis of damped nonrotating structures, 

The utility of the direct analysis of the state equation 

may depend on the computational efficiency of the inver- 

sion of <I> (after truncation), It should bo noted that is 

not orthogonal, i,e„ <i» »?M> 7 ', Foss (Ref, 36) has shown 
that when O' is zero, D ' is nonzero, and eigenvalues are 
distinct, certain weighted orthogonality relationships can 
be established, Foss notes the equivalence of Kqs, (157) 
and (140) with the first-order equation 


or, with Eq, (172) and truncation, 


A. 



0 


0 



(173) 
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This equation combines with the vehicle equation (169) 
and parallel equations for rotor and damper to obtain a 
complete simulation. Equation (173) is not applicable to 
the differential equations of an appendage on a rotating 
base (Eqs, 140 and 143 with G'-f O), When the base is 
not rotating, and specific damping forces are necessary 
(perhaps because of discrete dashpots in the system), the 
method of Foss is computationally superior, but mathe- 
matically equivalent, to the more general method used 
to derive Eq. (168), Both of these methods have two 
major undesirable features; 

(1) Except in the very special case for which G ' ' 0 
and D' is a polynomial in the symmetric matrices 
M' and K' (Ref, 30), the eigenvectors forming the 
columns of the transformation matrix 4* are com- 
plex, and thus the new* deformation variables Y 
introduced in Eq, (159) are also complex, Although 
physical interpretation of the complex coordinates 
(as conjugate pairs) has been provided, this inter- 
pretation is less immediate and probably less com- 
fortable for the engineer than would be a system of 
real coordinates, New computer programs might be 
necessitated by the simulation in terms of com- 
plex coordinates, whereas existing programs might 
suffice if a useful transformation to real coordinates 
can be found. 

(2) The final equation of motion of the appendage 
(Eq, 168 or Eq. 173) is a first-order matrix equation, 
Although controls engineers are finding increas- 
ing use for first-order equations, the structural 
dynamieist is more accustomed to working with 
second-order equations, Computation with existing 
computer programs may be precluded and inter- 
pretation of results may be impeded by the adop- 
tion of first-order equations. 


The general solution for Eq, (143) with A' ” 0 has 
already been found, since q is identified in Eq. (144) as 
the upper half of Q, and the solution for Q has been re- 
corded in real terms in Eq. (165). As noted previously, 
or in is zero for all m for undamped systems. If now the 
upper halves of 4' and r are denoted by the corresponding 
lower-case letters (as in Eq, 148 for <t»), the substitutions 


■vJ/W 



pm -5 


y m 1 

y m Am 


/ 


(17<i) 


permit q to be written from Eq, (165) in the form 


«n 


q " 2 { Cm (f" cos <r m ( - y m sin c T„t) 


III I 


D m (^'" sin (r m t 4- y m cos crmO } (175) 


The derivative q may be obtained by differentiating q or 
by inspecting the bottom half of Eq, (165), as 


mi 


() s { - Cmffm Sin <r m i + y m COS (i m t) 


m I 


+ Dm <7», (f n COS <T m t — Y w sin (TiiiO) (176) 


If now 7< h is defined as the 12n by 1 matrix 


4 , Modal analysis of undamped, rotating structures, 
When in Eq. (143) the acceleration matrix A' and the 
damping matrix D' are zero, an alternative to the trans- 
formation of Eq, (159) can he devised that results ulti- 
mately in a system of uncoupled, real, second-order scalar 
equations of motion, This method is inapplicable when 
damping is included in the mathematical model of the 
appendage structure, although the introduction of modal 
damping into the transformed equations is not precluded, 
The transformation to be described will uncouple the 
homogeneous second-order appendage equations even 
when the appendage is on a rotating base (so that G' / 0 
in Eq, 143), but this will not ho successful if a discrete 
dashpot is incorporated into the appendage or the ap- 
pendage support structure, since then D' 6 in Eq, (143), 


C, cos <7,1 + D x sin <7 { t 

Cy COS <7y t I D« sin <7,4 




Cm, cos <7,1,4 T Do,, sin <7<w4 

sin a,# 4- P, cos <7,1 
< 

! 

1 

C u „ sin <7(o4 4 Don cos <7 ui4 


(177) 


then by inspection of Eqs, (165) and (177), one may write 

Q « m (178) 


JPt TECHNICAl REPORT 33- 1339 


47 



whore P is the transformation matrix 


liU ' 

y 1 

{jj~ j , • • j I/, 0 ** 

, p „ - - - 

r> 

u 1 

! 

1 

1 

! 1 

(179) 

_ — cr t y' 

— (T..,y" 1 ' * ’ 1 o r ,i„Y (Vn 


i 

' ' l Vi\nV H 



By inspection of the matrix Z h and its time derivative, one 
may wife 


in the real variables s$i, ‘ * ' ,So„, which is equivalent to 
Eq, (1-13). In matrix terms, this result is 


0 


I <n 0 
I \ 


Z 1 ' 


\ 

T 


0 




<*1 


0 


0 


Van 


7l { \ 


0 j (r 

(T | 0 


z h 


(.180) 


where or has been defined as the On by On diagonal matrix 
of natural frequencies <q, * * * ,<r UH , 


-H irs 


0 


(183) 


where 


<r'\ 0 

0 cfili 


To accomplish a direct transformation from q in 
Eq, (143) to 3 in Eq, (183), rewrite the transformation in 
Eq, (1 44) as 

<I^W\0]Q (184) 


The first On elements of Z h may be assigned the sym- 
bols • • * , Sun, so that Eq, (177) may be written as 


where the identity matrix E and the null matrix 0 are of 
dimension Go by On, Substitution of Eqs, (178) and (181) 
into Eq, (184) produces 


■v. *n 

**\ 
t 



* 

t 

Sum 


r * 1 


r-— ; 

~~ w — 

£uM 


iT 1 

(* A* 


(181) 


q - [E\ 0] P 



or, from Eq, (179), 



q ~ + yor 



(180) 


which defines the On by 1 column matrix % and the On by 
On diagonal matrix or 1 whose elements are the frequency 
reciprocals I/m, • • • , l/<r UM , In terms of Si, • * 1 ,s««» 
Eq, (180) becomes the set of scalar equations 

i « 

<•»* rr: 

<\'\ — <V1 

t 

_ ' 

*"*• <vUM 

'i\ = 

t 
\ 

U Q ^ 

o'an “ <*oh*uh 

The first set of On scalar equations in Eq, (188) is trivial, 
but the second set is a group of On uncoupled equations 



where and y are the On by On matrices 

^ ~a 1 > * 

and 

yES^y* 1 ’ ' y an ] 

A transformation of the character of Eq, (185), which 
relates old coordinates to a combination of new coordi- 
nates and their derivatives (or, more often, their gen- 
eralized momenta), is known as a contact transformation. 
In contrast, a transformation that relates one set of gen- 
eralized coordinates directly to another is called a point 
transformation , As noted in Whittaker (Ref, 24, p, 195, 
footnote), there exists no point transformation that takes 
equations of the structure of Eq, (143) with 
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into uncoupled second-order equations such as Eq, (183), 
unless G' is also zero, The e.dstenee of a contact trans- 
formation that accomplishes this objective when G' /- 0 
is noted in Whittaker (p. -127), where equations of motion 
are expressed in terms of Hamilton’s canonical equations, 


Only when L is zero is the solution for Z that given as Z h 
by Eq. (177), 

When P is nonsingular (a condition not always assured), 
Eq, (188) may be written as 


Although Eq, (1-13) has been the starting point in the 
quest for useful transformations, it is the inhomogeneous 
counterpart, Eq, (1-10), that must actually be employed 
in the simulation, In terms of first-order equations, it is 
Eq, (157) and not Eq, (145) that must be used. Direct 
substitution of the transformation in Eq, (186) into the 
second-order inhomogeneous Eq, (140) does not appear 
to he fruitful. Instead the transformation 

Q ~ FZ (187) 


Z = P X BPZ + P l L (189) 


The matrix P X BP has been shown in the vontext of the 
analysis of the homogeneous equation (see Eq. 180) to 
have the value 


P 1 BP 


0 t" 

”"""T 

” (T i 0 


(190) 


is used in the first-order inhomogeneous Eq, (157), where 
P is a matrix of constants available from the eigenvalues 
and eigenvectors of B as in Eq, (179), and Z is a new 
matrix of unknowns. The result is, from Eqs, (157) and 
(187), 

PZ BPZ ( L (18S) 


and this product of constant matrices cannot change its 
value when the equation becomes inhomogeneous. 

The inverse of P as required in Eq, (189) is expressed 
and computed most easily in terms of its partitioned ele- 
ments, according to the formula (Ref, 37, p. 640) 


P ' 



'PdP.t (P 


PnP^Pyj_ 

u “ P^P,lP, t )-‘ 


! (p k - p«r„'P «) 1 . 


(191) 


where 



' Pn ! P« " 


t | y 

Pa» 

l 


~ 1 


L p.u I p a , J 


_ ~ vy 1 (rp _ 


With the identity typified by 



(193) 


the inverse of P simplifies to 


(^ + y\j/~ x y)~ x | — \}i~ x y ($ -f- x y) x <j l 
' y (p 4 y$ x y) x I (^ + y\f> x y) x a 1 


(194) 
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In expanded form, Eq. (189) may now be written (see Eq. 157 for L expansion) in the form 


" “■ 


1 cr t 0 


“ — 

z t 


o 1 \ 


lx 

* 


1 0 (Jri i 



• 


— cr t 0 | 



Zu„ 


\ ! 0 


z l2M 

- 


- 0 — cron 1 


J 


(^ + y^' 1 y)-l 

— if/" 1 y (^ 4- yij/~ x y)~ l <r~ l ~j 

* 0 

_ i’ 1 y 4" r f * yY x 

0 4- yf 1 y)- 1 or' 1 J 

L'_ 


(195) 


In the analysis of the homogeneous equation, with V 
absent, the elements of the column matrix Z, labeled Z\ 
can be identified as in Eq, (181). This permits the scalar 
equations obtained from the upper half of Z to be dis- 
carded as trivial identities (see Eq. 182), and leaves 
the lower half of Z to provide nontrivial, uncoupled 
second-order equations in the new variables 
(Eq. 183). This identification and resolution into second- 
order equations does not appear to be possible in the in- 
homogeneous case (Eq. 195), unless — f *y (i> 4 #'V)"V 1 
is zero. Note that the inverted matrix in parentheses can- 
not be zero, or every term in P 1 is zero, and the applied 
forces and torques in V have no dynamic consequences. 
The most useful condition under which Eq. (195) permits 
this identification of the elements of Z as in Eq, (181) is 
the condition y ~ 0. This corresponds to the special case 
of real eigenvectors, which implies the absence of O' in 
Eqs. (1-13) and (140). This is then the classic problem of 
the vibration of a structure about a state of inertial rest, 
to which attention will return shortly. 

In general, the 12 n first-order equations in Eq. (195) do 
not admit of reduction to second-order equations in z as 
in the homogeneous case (Eq. 183). This reduction does 
serve, nonetheless, to motivate the coordinate truncation 
of the deformation variable matrix Z in the inhomogene- 
ous Eq. (195). 


z>n C'hi cos cr m t T D n i sin <T„i t 

z m . O’,,) Cm sin o’mt 4" <r ,1) Dm COS (Tmt 

then substitution into Eq. (186) yields the solution for q 
in Eq. (175), for which physical interpretation has been 
provided. 

The introduction of energy-dissipation capability into 
the mathematical model can also be accomplished in the 
homogeneous case after transformation to the modal- 
coordinate matrix 3, Equation (183) is simply replaced by 

3 4- 2 4 cr5 = 0 (197) 

where the overbars indicate truncation to N coordinates, 
and 


Ci 0 



In the inhomogeneous case, which is of primary inter- 
est, the matrix Z in Eq. (195) can be written in partitioned 
terms as 



With the use of frequency criteria, the truncation of 
the 6u by 1 matrix % in Eq, (183) is completely straight- 
forward; one simply ignores all but ' ' ' ,s.v of the 
elements in s when frequencies tr u ' ’ 1 , vx are of inter- 
est. The coordinates %, * * * ,s.y are modal coordinates 
in the broad sense in which that term is used in this 
report. If the solutions to the scalar equations in Eq, (1S3) 
are written as 


Z l 


z= 



(198) 


where Z l and 2- are 6n by 1 matrices. The transformation 
derived as Eq. (1S6) for the homogeneous case then 
becomes 


q - ^Z 1 4- y Z- (199) 
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